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The  use  of  antenna  arrays  is  a promising  technology  for  improving  performance  of 
all  types  of  cellular  networks.  This  dissertation  investigates  the  use  of  multiple  antenna 
sensors  for  parameter  estimation  and  detection  in  Direct-Sequence  Code-Division 
Multiple- Access  (DS-CDMA)  communication  systems.  Our  focus  is  on  the  timing 
estimation  that  is  essential  for  any  type  of  DS-CDMA  system.  In  order  to  increase  the 
performance  and  capacity  of  the  DS-CDMA  system,  the  use  of  an  adaptive  antenna  array 
and  diversity  combiner  has  received  much  attention,  but  there  has  been  less  attention  to 
the  code  timing  estimation  problem  when  using  multiple  antenna  sensor  arrays. 

We  describe  a general  spatial  channel  model  for  exploiting  multiple  antenna  sensors 
at  the  base  station.  Tap  weight  vector  estimator  for  beamforming  or  diversity  combining 


IX 


is  proposed  and  evaluated.  We  derive  the  Cramer-Rao  Bound  (CRB)  for  the  parameter 
estimation  problem  when  using  multiple  antenna  sensors.  The  code  timing  estimation 
problem  with  multiple  antenna  sensors  is  introduced.  The  adaptive  antenna  array  and 
diversity  combiner  based  code  timing  estimation  schemes  are  proposed  and  studied.  The 
proposed  approach  for  estimating  the  code  timing  is  shown  to  improve  the  performance 
of  a single  antenna-based  timing  estimator  in  both  fading  and  nonfading  conditions  when 
we  use  an  existing  CDMA  system’s  standard  power  control  technique.  The  problem  of 
estimating  multipath  timings  of  DS-CDMA  signals  when  using  a multiple  antenna  sensor 
array  which  exploits  spatial  diversity  is  also  considered  on  a frequency  selective  fading 
channel 

The  minimum  mean-squared  error  (MMSE)  receiver  is  a very  attractive  receiver  for 
the  DS-CDMA  signal  detection.  We  study  the  effect  of  the  multiple  antenna  sensors  on 
the  MMSE  detector  over  fading  channels. 


CHAPTER  1 

INTRODUCTION  AND  BACKGROUND 


In  this  chapter,  the  background  of  direct-sequence  code-division  multiple-access 
(DS-CDMA)  and  multiuser  detection  will  be  presented.  This  chapter  will  introduce 
timing  estimation  based  on  a single  antenna  and  multiple  antenna  sensors.  An  outline  of 
the  dissertation  will  conclude  this  introductory  chapter. 

1.1  Direct-Sequence  Code-Division  Multiple- Access  (DS-CDMA) 

Spread  spectrum  communications  is  a digital  communication  technique  whereby  the 
transmitted  signal  occupies  a much  larger  bandwidth  than  the  minimum  bandwidth 
required  to  transmit  the  signal.  This  bandwidth  spreading  is  accomplished  by  means  of 
code  sequences  known  to  both  the  transmitter  and  receiver.  Spread  spectrum 
communications  was  developed  in  military  applications  for  its  robustness  to  a wide 
variety  of  interference.  There  are  two  major  spread  spectrum  techniques.  One  is  Direct- 
Sequence  Spread- Spectrum  (DS-SS)  where  the  transmitted  signal  is  modulated  by  a high- 
rate  pseudorandom  sequence  known  as  a spreading  sequence  and  the  other  is  Frequency- 
Hop  Spread-Spectrum  where  the  carrier  frequencies  of  the  individual  users  are  assigned 
in  a pseudorandom  pattern. 

Since  spread  spectrum  signals  are  inherently  resistant  to  interference,  its  system  can 
provide  multiple  access  capability  for  commercial  communication.  This  is  referred  to  as 
Code-Division  Multiple- Access  (CDMA).  We  shall  concentrate  on  direct-sequence 
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CDMA  (DS-CDMA).  The  DS-CDMA  transmitter  multiplies  each  user’s  signal  by  a 
unique  spreading  sequence.  The  signal  sent  by  a desired  user  can  then  be  demodulated  by 
correlating  the  received  signal  with  the  code  sequence  of  the  desired  user.  Only  the  one 
signal  using  the  user’s  code  sequence  is  despread  properly.  This  demodulator  that  is 
widely  used  in  practice  is  referred  to  as  the  conventional  single-user  detector.  DS-CDMA 
has  received  much  attention  for  its  application  to  multi-user  systems  such  as  mobile 
cellular  and  personal  communication  systems. 

A U.S.  digital  cellular  system  based  on  DS-CDMA  has  been  standardized  as 
Interim  Standard  95  (IS-95)  by  the  U.S.  Telecommunications  Industry  Association.  The 
IS-95  standard  for  CDMA  cellular  systems  was  published  in  1992. 

1.2  Multiuser  Detection 

Conventional  detection  of  DS-CDMA  signals  is  performed  by  matched  filtering  the 
received  signal  against  the  desired  user’s  code  waveform.  Due  to  the  lack  of 
orthogonality  between  user’s  code  waveforms,  a portion  of  the  extra  users’  signal  power 
appears  to  the  receiver  like  interference  after  matched  filtering.  That  is,  in  a conventional 
CDMA  system,  all  users  interfere  with  all  other  users.  This  multiuser  interference  is 
referred  to  as  Multiple  Access  Interference  (MAI).  As  the  number  of  interfering  users 
increases,  the  error  performance  will  be  degraded.  Even  if  the  number  of  users  is  not  too 
large,  the  strong  nearby  interfering  users  may  swamp  a weak  desired  user’s  signal  far 
away  from  the  receiver.  This  is  known  as  the  near-far  effect.  One  technique  to  solve  the 
near-far  problem  is  to  use  transmit  power  control  such  as  is  used  in  the  IS-95  CDMA 
system.  That  is,  each  mobile  continually  adjusts  its  power  so  that  the  received  power  at 
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the  base  station  from  all  users  is  equal.  Among  the  benefits  of  the  use  of  power  control, 
the  following  can  be  cited.  In  a mobile  urban  environment  with  typical  cell  sizes, 
imbalances  of  the  received  powers  of  all  users  can  be  as  high  as  90  dB.  Power  control 
solves  the  near-far  problem  for  such  an  enormous  dynamic  range.  By  expending  no  more 
power  than  that  required  by  the  receiver,  it  can  save  the  power  consumption  of  portable 
transmitters.  So,  this  extends  the  battery  life  of  portable  transmitters.  It  can  also,  to  some 
extent,  alleviate  the  effects  of  fading. 

The  conventional  single-user  detector  is  optimum  for  a single  user  when  the 
received  signal  is  corrupted  by  additive  white  Gaussian  noise  alone.  It  is  not  however  the 
optimum  detector  under  the  existence  of  MAI.  Because  of  the  MAI,  a better  detection 
strategy  that  is  called  multiuser  detection  is  needed.  The  optimum  multiuser  detector  [3] 
or  the  maximum  likelihood  sequence  detector  uses  information  about  multiple  users  to 
better  detect  each  individual  user  and  is  shown  to  consist  of  a bank  of  single  user  matched 
filters  followed  by  a Viterbi  algorithm.  Unfortunately,  this  detector  is  much  too  complex 
for  practical  DS-CDMA  systems.  It  also  requires  a great  deal  of  side  information  like  the 
number  of  users,  the  code  waveforms  of  all  users,  the  received  amplitudes  of  all  users  and 
the  code  timings  of  all  users.  Therefore,  several  suboptimal  multiuser  detectors  [4-12] 
that  are  more  feasible  to  implement  have  been  proposed.  These  multiuser  detection 
schemes  are  resistant  to  the  near- far  problem  and  have  greatly  reduced  complexity. 
Among  the  near-far  resistant  suboptimum  receivers,  we  can  cite  the  decorrelating 
receiver  and  the  minimum  mean  square  error  (MMSE)  receiver.  The  decorrelating 
detector  [4,5]  applies  the  inverse  of  the  cross-correlation  matrix  of  the  matched  filter 
outputs  to  the  conventional  detector  output  in  order  to  decouple  the  data.  It  requires  less 
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side  information  like  the  code  waveforms  of  all  users  and  the  code  timings  of  all  users 
than  that  of  the  optimum  detector.  The  MMSE  receiver  [9-12]  has  received  significant 
attention  for  use  in  detection  of  DS-CDMA  signals.  It  is  a single  user  detector  in  that  it 
only  demodulates  the  bit  stream  for  one  user.  It  uses  an  adaptive  algorithm  to  obtain  the 
optimum  solution  for  its  coefficients.  It  requires  little  side  information.  Only  code  timing 
of  the  desired  user  must  be  acquired. 

1.3  Timing  Estimation 

Multiuser  receivers  assume  the  knowledge  of  some  parameters,  including  the 
propagation  delay  for  each  user,  received  powers  and  carrier  phases.  In  this  work  we  will 
focus  on  the  timing  estimation  problem  in  DS-CDMA  communication  systems  because  it 
is  difficult  to  obtain  accurate  propagation  delay  estimates.  It  should  be  emphasized  that 
the  chip  tracking  is  essential  for  any  type  of  DS-CDMA  system  and  that  tight  code 
synchronization  within  a small  fraction  of  one  chip  duration  is  required  for  reliable 
operation  of  the  conventional  detector.  The  capacity  of  a DS-CDMA  communication 
system  employing  the  matched  filter  approach  in  the  presence  of  MAI  was  shown  to  be 
limited  by  the  timing  estimation  problem  [14]. 

The  question  is  how  much  tighter  the  propagation  delay  estimation  error  tolerable 
for  the  conventional  detector  must  be  for  multiuser  detection.  In  [15,16],  the  analysis  of 
the  sensitivity  of  a near- far  resistant  receiver  to  propagation  delay  estimation  error  has 
been  addressed.  The  decorrelating  receiver  is  shown  to  lose  its  near-far  resistance  in  the 
presence  of  imperfect  propagation  delay  estimates. 
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The  problem  of  obtaining  more  accurate  and  efficient  timing  has  received  more 
attention  recently.  In  [18,19,23],  the  minimum  mean  squared  error  (MMSE)  receiver 
based  propagation  delay  estimator  is  presented.  Schemes  for  the  joint  acquisition  and 
demodulation  of  CDMA  signals  are  presented  in  [21].  Timing  estimation  using  the 
MUSIC  algorithm  was  discussed  in  [20]  and  [22].  The  MUSIC  based  timing  estimator 
requires  more  computational  complexity  due  to  the  eigendecomposition  operations  than 
that  required  by  the  correlator  type  estimators.  In  [17],  an  efficient  timing  estimator  using 
a known  training  sequence  was  developed.  The  timing  estimators  for  code  acquisition 
mentioned  above  are  based  on  a single  antenna  element. 

1.3.1  Single  Antenna-based  Timing  Estimator 

1.3.1. 1 Sliding  correlator  T21 

This  is  the  conventional  approach  to  code  acquisition  whereby  the  received  signal  is 
correlated  with  time  delayed  versions  of  the  known  code  sequence  and  the  timing 
estimate  is  simply  the  value  of  the  time  delay  which  maximizes  the  correlation.  The 
sliding  correlator  is  computationally  simple  and  it  is  optimal  for  a single  user  in  the 
presence  of  white  Gaussian  noise  only,  but  can  be  highly  suboptimal  in  the  presence  of 
MAI. 

1.3. 1.2  MMSE-receiver-based  timing  estimator 

This  approach  was  developed  by  Smith  and  Miller  [18,19,23]  and  is  based  on  a 
MMSE  detector.  The  adaptive  filter  tap  weights  are  chosen  to  minimize  the  mean  squared 
error  between  the  training  sequence  and  the  filter  output.  Once  the  MMSE  receiver  tap 
weights  are  computed,  the  timing  can  be  estimated  from  these  weights  using  a correlator 
type  approach. 
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1.3. 1.3  MUSIC  timing  estimator  f201 

The  space  of  all  received  vectors  can  be  broken  into  a signal  space  and  a noise 
space  which  can  be  found  from  the  eigendecomposition  of  the  autocorrelation  matrix  of 
the  received  vector.  Since  the  known  code  sequence  is  contained  in  the  signal  subspace,  it 
must  be  orthogonal  to  the  noise  subspace.  Thus,  we  can  find  the  timing  delay  of  the 
desired  user  by  finding  the  value  of  timing  delay  for  which  the  known  code  sequence  is 
orthogonal  to  an  estimate  of  the  noise  subspace. 

1.3. 1.4  LSML  timing  estimator 

A large  sample  maximum  likelihood  (LSML)  timing  estimator  [17]  models  the 
known  training  sequence  as  the  desired  signal  and  all  other  signals  including  the 
interfering  signals  and  additive  noise  as  unknown  colored  Gaussian  noise  which  is 
uncorrelated  with  the  desired  signal.  The  problem  then  becomes  one  of  estimating  the 
timing  of  a single  desired  signal  in  noise  with  unknown  autocorrelation.  The  LSML 
timing  estimator  will  be  given  in  Chapter  5. 

1.3.2  Adaptive  Antenna  Array  and  Diversity  Combiner 

The  goal  of  wireless  communication  systems  is  to  provide  a wide  variety  of 
wireless  communication  services  via  the  efficient  use  of  limited  available  frequency 
spectrum  for  many  users.  More  and  more  users  are  requesting  those  services.  The  current 
design  trend  for  those  systems  is  such  that  the  coverage  area  is  divided  into  smaller 
regions  or  cells.  A large  number  of  small  cells,  each  served  by  a base  station,  result  in 
increased  base  station  equipment,  increased  networking  and  increased  handoff.  When  a 
mobile  user  crosses  the  boundary  between  two  cells,  its  communication  channel  is 
switched  to  the  base  station  in  the  new  cell.  This  is  known  as  handoff.  For  a system 
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capacity  increase,  IS-95  CDMA  system  uses  120°  sectorization  by  using  directional 
antennas  at  each  base  station. 

The  use  of  spatial  filtering  (beamforming)  or  diversity  combining  approaches  via 
multiple  antenna  sensor  arrays  is  considered  to  be  a promising  technology  for  providing 
an  increased  capacity  and  performance  for  cellular  mobile  and  personal  communication 
services  [28-40].  The  adaptive  antenna  array  (or  adaptive  beamformer)  exploits  the 
available  information  such  as  spatial  dimension  or  reference  signal  in  signal  processing 
for  wireless  communication  networks.  Its  beams  in  the  antenna  pattern  are  artificially 
formed  in  the  direction  of  the  desired  signal  while  minimizing  the  influence  of  the  others. 
Its  beams  are  formed  by  a linear  combination  of  the  array  output  that  is  a beamformer 
output.  It  exploits  direction-dependent  phase  shifts  between  any  adjacent  antenna 
elements.  The  other  technique  used  to  improve  performance  over  a fading  channel  is  to 
provide  the  receiver  with  multiple  versions  of  the  same  transmitted  signal  each  with 
statistically  independent  fading.  One  of  those  diversities  can  be  achieved  through  space 
by  using  several  receiving  antennas,  which  is  called  spatial  diversity.  A spatial  diversity 
combiner  combines  the  received  signals  from  multiple  antenna  sensors  sufficiently 
separated  in  space  with  uncorrelated  fading.  From  now  on,  we  will  omit  the  word  spatial 
from  the  term  ‘spatial  diversity  combiner.’ 

The  diversity  combiner  is  different  from  the  adaptive  antenna  array  in  the  sense  that 
it  does  not  process  any  spatial  information  like  the  direction  of  arrival  (DO  A)  of  the 
received  signal.  It  combines  the  received  signals  to  increase  the  signal  level  without 
affecting  the  individual  antenna  pattern.  On  the  other  hand,  the  adaptive  beamformer 
exploits  the  differential  phase  between  different  antennas  to  change  the  antenna  pattern  of 
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the  whole  array.  After  it  combines  the  received  signals,  the  whole  of  the  array  has  a 
single  antenna  pattern.  In  Chapter  2,  we  will  show  that  as  the  number  of  DOAs  of  the 
incoming  signals  at  the  base  station  goes  to  infinity,  the  system  model  for  the  adaptive 
antenna  array  becomes  that  of  a diversity  combiner.  The  key  to  the  adaptive  beamforming 
and  diversity  combining  techniques  is  the  determination  of  the  combining  weight  vector. 

The  main  advantages  of  using  the  adaptive  antenna  array  are  as  follows.  The  first 
one  is  to  improve  the  signal  to  interference  ratio,  thus  reducing  the  transmit  power 
required  by  the  mobile  handset  for  the  reverse  link,  Le.,  mobile  to  base  station.  The 
second  one  is  to  increase  the  number  of  active  users  for  a given  bit  error  rate  (BER) 
quality  threshold.  The  third  one  is  to  increase  the  range  of  the  base  station  receiver  and 
thus  the  cell  size.  The  fourth  one  is  to  permit  a less  stringent  form  of  reverse  link  power 
control  while  maintaining  acceptable  BER  performance.  Diversity  combining  is  used  to 
mitigate  the  effects  of  fading  of  the  desired  signal  and  suppress  the  interfering  signals  and 
thereby  increase  output  signal  to  interference  plus  noise  ratio  (SINR). 

In  order  to  increase  the  performance  and  capacity  of  the  DS-CDMA  system,  the  use 
of  an  adaptive  antenna  array  and  diversity  combiner  has  received  much  attention,  but 
there  has  been  much  less  attention  to  the  code  timing  estimation  problem  when  using 
multiple  antenna  sensors  array  for  DS-CDMA  communication  systems. 

1.3.3  Previous  Work  on  Code  Acquisition  via  Multiple  Antenna  Array 

Previous  work  which  has  dealt  with  the  acquisition  problem  of  spread  spectrum 
(SS)  signals  using  multiple  antenna  arrays  includes  [24],  [25],  [26],  and  [27].  In  [24], 
Dlugos  and  Scholtz  consider  an  array  for  acquisition  of  SS  signals  in  a single-user 
scenario  using  the  maximum  likelihood  (ML)  estimator,  but  these  techniques  cannot  be 
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used  when  a large  number  of  users  are  present.  In  [25],  the  MUSIC  algorithm  that 
exploits  spatial  diversity  at  the  receiver  using  antenna  arrays  to  estimate  the  code  timing 
is  discussed.  In  [26]  and  [27],  the  LSML  code  timing  estimator  in  [17]  is  extended  to  the 
case  of  using  multiple  antenna  sensor  array  which  exploits  receiver  diversity.  But,  the 
techniques  described  in  [24]-[27]  are  based  on  a diversity  model  and  hence  have  no 
beamforming  process  to  suppress  the  interfering  signals  to  help  in  estimating  code 
timing. 

1.3.3.1  MASE  algorithm 

A multiple  antenna  sensors  based  estimator  (MASE)  [26]  is  an  extension  of  the 
LSML  code  timing  estimator  in  [17]  which  is  based  on  a single  antenna  to  a receiver 
diversity  DS-CDMA  system  employing  multiple  antenna  arrays.  The  MASE  algorithm  is 
an  asymptotic  (for  large  data  samples)  maximum  likelihood  estimator  that  is  derived  by 
modeling  the  known  training  sequence  as  the  desired  signal  and  all  undesired  signals 
including  the  MAI  and  the  additive  noise  as  unknown  colored  Gaussian  noise.  The 
authors  of  [26],  using  the  MASE  algorithm,  consider  only  the  case  of  a single  path  flat 
fading  channel  and  do  not  consider  the  case  of  frequency  selective  fading  due  to 
multipath  propagation.  Moreover,  the  amount  of  computations  required  by  MASE  is 
much  more  than  that  required  by  the  LSML  timing  estimator.  The  MASE  algorithm  will 
be  discussed  further  in  Chapter  5. 
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1.4  Dissertation  outline 

This  dissertation  investigates  how  multiple  antennas  can  be  used  for  parameter 
estimation  and  detection  in  DS-CDMA  communication  systems.  Timing  estimation  and 
MMSE  detection  based  on  multiple  antennas  which  can  discriminate  between  spatially 
separated  users  and  provide  spatial  diversity  are  proposed.  Chapter  2 describes  a general 
spatial  channel  model  including  spatial  information  such  as  DO  As  and  spatial  diversity. 
In  Chapter  3,  two  techniques  for  computing  the  tap  weights  of  the  adaptive  antenna  array 
or  diversity  combiner  are  discussed.  An  adaptive  antenna  array  or  diversity  combiner  can 
suppress  the  interfering  signals,  Le.,  MAI  and  hence  the  output  signals  of  an  adaptive 
antenna  array  or  the  diversity  combiner  have  less  interference.  Since  the  timing 
acquisition  is  affected  by  the  multiple  access  interference  [14],  more  accurate  and  MAI- 
resistant  timing  estimates  could  be  obtained  if  we  use  the  output  data  of  the  adaptive 
antenna  array  and  diversity  combiner.  The  Cramer-Rao  Bound  (CRB)  is  derived  for  the 
parameter  estimation  problem  based  on  multiple  antennas  over  an  AWGN  channel  and 
fading  channels  in  Chapter  4.  In  Chapter  5,  we  apply  the  LSML  timing  estimator  [17]  to 
the  output  data  of  the  adaptive  antenna  array  or  diversity  combiner  to  estimate  the 
propagation  time  delay  of  the  desired  user  in  multi-user  scenarios.  It  is  shown  that  the 
adaptive  antenna  sensor  array  and  diversity  combiner  make  the  LSML  estimator  achieve 
improved  performance  in  terms  of  the  root  mean  squared  error  (RMSE)  and  probability 
of  correct  acquisition  on  an  AWGN  channel  and  time-varying  fading  channels, 
respectively.  This  means  that  it  offers  better  performance  when  the  number  of  users  is 
large.  We  extend  the  MASE  algorithm  to  a frequency  selective  fading  channel  for  joint 
estimation  of  multipath  timings  of  a desired  user’s  signal.  The  acquisition  performance  of 
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the  MUSIC  timing  estimator  with  a single  antenna  is  not  good  for  a small  number  of 
observation  bits  [47]  and  also  for  low  signal  to  noise  ratios  (SNR)  [22].  By  exploiting 
spatial  diversity  via  multiple  antennas  its  acquisition  performance  even  for  a small 
number  of  observation  bits  or  low  SNRs  could  be  improved  as  compared  to  that  of 
MUSIC  with  a single  antenna.  So,  the  MUSIC  timing  estimator  in  [47]  is  extended  to  the 
case  of  using  a multiple  antenna  array  for  estimating  multipath  timings  in  frequency 
selective  fading  channels.  The  multiple  MMSE  detectors  with  multiple  antenna  sensors 
are  introduced  in  Chapter  6.  We  describe  how  the  spatial  diversity  via  a multiple  antenna 
array  will  help  the  performance  of  the  MMSE  detection.  Chapter  7 gives  a summary  of 
the  results  presented  and  some  directions  for  future  research. 


CHAPTER  2 
SYSTEM  MODEL 


In  this  chapter  we  will  develop  a general  channel  model  that  covers  all  special  cases 
in  which  this  dissertation  is  interested.  In  the  work  throughout  this  dissertation,  Qw , (•)* , 

and  (*)r  denote  the  complex  conjugate  transpose,  complex  conjugate,  and  transpose, 
respectively. 


2.1  Continuous-Time  System  Model 

We  consider  an  asynchronous  K-user  DS-CDMA  system.  The  modulation  scheme 
is  BPSK  with  bit  duration  T and  chip  duration  Tc  - T/N , where  N is  an  integer.  The  code 

waveforms  have  unit  amplitude  and  are  assumed  to  be  rectangular  and  periodic  with 
period  T.  Asa  general  rule,  a subscript  k implies  that  the  quantity  is  due  to  the  kth  user. 
For  instance,  a period  of  the  Jtth  user’s  code  waveform  is  denoted  by  bk  (t) , where 

N 

bk  ( t ) = 0 for  t€  [0, T ) . That  is,  bk  ( t ) = '£Jck  (rc)nr  ( t-{n  - 1 )TC),  where 

n=\ 

ck  (n)  e {+1,-1}  and  nr  (t)  denotes  a unit  rectangular  pulse  over  the  chip  period  [0, Tc) . 

The  baseband  signal,  sk  (t) , is  formed  by  pulse  amplitude  modulating  the  data 
stream,  dk  (m)  e {+1,-1} , with  a period  of  the  code  waveform,  Le., 

M 

sk(t)  = '£idk(m)bk(t-mT ).  (2.1) 

m= 0 
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The  transmitted  signal  is  formed  by  multiplying,  sk  (t) , with  the  carrier  ^2Pk  co s(wct) , 
where  Pk  is  the  transmitted  power.  We  assume,  without  loss  of  generality,  that  P1  = 1 
and  Tc  = 1 . 

In  a mobile  communication  system,  a signal  transmitted  into  a wireless  channel 
interacts  with  the  propagation  environment  with  physical  objects  in  a complicated 
process  that  involves  reflection,  diffraction  and  scattering.  There  exist  many  signal 
components  from  these  complex  interactions  and  this  situation  is  referred  to  as  multipath 
propagation.  The  effect  of  time-variant  multipath  characteristics  of  the  channel  can  cause 
the  variations  in  the  received  signals  amplitude,  phase  and  direction  of  arrival  (DOA). 
This  phenomenon  is  known  as  multipath  fading. 

In  Figure  2.1,  we  have  illustrated  a simplified  wireless  multipath  propagation 
scenario.  Since  each  signal  component  arrives  at  the  base  station  over  a different  path 
environment,  it  has  own  amplitude  ak  r,  carrier  phase  shift  Qk  r,  propagation  time  delay 

zk  r from  mobile  to  base  station  and  DOA  <pk  r of  the  rth  signal  component  of  the  kth 

mobile.  In  general,  each  of  these  signal  parameters  will  be  time  varying. 

The  classical  channel  model  proposed  by  Turin  [41]  for  the  kth  mobile  user  is  a 
time-varying  filter  with  impulse  response  hk  ( z,t ) 


where  8(t ) is  the  Dirac  delta  function,  and  Rk  (t)  is  the  number  of  paths.  The  amplitude 
of  the  fading  process  ak  r(t ) is  Rayleigh  distributed  and  the  phase  of  the  fading  process 
dk  r 0)  is  uniformly  distributed. 


(2.2) 
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(a)  side  view 


Figure  2. 1 : Multipath  propagation  channel  model 

The  channel  model  in  (2.2)  does  not  consider  the  DOA  of  each  multipath 
component  shown  in  Figure  2.1.  We  consider  the  narrowband  signals  that  refer  to  the  fact 
that  each  signal  has  a small  bandwidth  relative  to  the  carrier  frequency  wc . If  the 
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narrowband  signals  are  considered  to  be  impinging  at  an  array  of  P antenna  elements,  the 
Pxl  vector  spatial  channel  impulse  response  by  including  the  DOA  parameter  may  be 
represented  as 

RJJ) 

r= 1 


where  d>(<pt  ,)  is  the  array  response  vector  to  a signal  from  direction  q>k  r which  is  a 
function  of  the  array  geometry. 

To  simplify  the  discussion,  we  assume  that  the  mobile  is  far  afield  and  lies  in  the 
same  plane  as  the  antenna  array.  In  this  case,  the  propagation  time  delay  At jfj  (t)  from  a 

reference  point  to  sensor  p (p  = 1,  2, , P)  of  the  antenna  array  depends  only  on 

azimuth  (pkr  ( t ) . Then,  the  array  response  vector  is  given  by 


®{(Pk,r)  = 


(2.4) 


A uniform  linear  array  (ULA)  is  the  array  of  P identical  sensors  uniformly  spaced 
on  a line  as  illustrated  in  Figure  2.2.  By  taking  the  first  antenna  sensor  as  a reference 
point,  the  array  response  vector  is  given  by 


<D(<pt  )=[l  ^2(i»-i x«< /*)««*., (*)]r  ^ 5' 

where  d is  the  distance  between  two  consecutive  antenna  sensors  and  X is  the  signal 
wavelength.  In  this  case,  we  can  define  the  parameter  lik  r{t)  = 2(jtd / A) sin  (pkr{t)  that  is 


the  direction-dependent  phase-shift  between  any  two  adjacent  array  elements  for  a ULA 
where  (pk  r is  uniformly  distributed  in  [~n/2,n/2)  . 

A change  in  the  array  geometry  would  affect  only  the  functional  form  of  the  relative 
direction-dependent  phase-shift  between  the  array  sensors.  If  we  could  consider  the  array 
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consisted  of  identical  and  omnidirectional  elements  which  are  uniformly  distributed  over 
the  circumference  of  a circle  of  radius  r in  the  xy  plane,  as  shown  in  Figure  2.3,  the  phase 

.2 nr  ( M 2k\ 

J oo«  f ( t ) 

difference  between  any  two  adjacent  sensors  is  given  by  e 1 ' p > [66]  where 

azimuth  angle  q>k  is  measured  counterclockwise  from  the  x axis.  That  is,  the  array 
response  vector  is  given  by 


0(<P*.r)  = 


r(O) 


.2  nr  f 

/Tci*" 


(D- 


(P~ D2* 


(2.6) 


Unlike  the  ULA,  the  uniform  circular  array  (UCA)  can  resolve  incident  signals  over  the 
DO  A interval  [0,  lie). 


Figure  2.2  : The  uniform  linear  array  scenario 
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Figure  2.3  : Uniform  Circular  Array  Geometry 

In  this  work,  we  assume  that  the  number  of  paths  and  their  propagation  delays  and 
DOAs  are  slowly  varying  in  time  compared  to  the  observation  time  for  the  received 
signal  In  fact,  the  number  of  paths  and  their  propagation  delays  and  DOAs  are 
considered  to  be  fixed,  Le.,  Rk  ( r ) = Rk , rkr(t)  =rk<r  and  (pkr(t)  = <pk  l.  for  re  [0,  MT). 

We  also  assume  that  clocks  of  the  transmitter  and  receiver  for  one  path  of  the  kth  user’s 
signal  have  been  aligned  to  roughly  within  a bit  interval  and  the  delay  spread  of  the 
channel  is  within  several  chip  intervals.  Hence,  we  consider  only  the  relative  propagation 
delay,  that  is,  rk  r e [0,r).  Furthermore,  the  fading  processes  are  assumed  to  be  wide- 
sense  stationary  and  to  vary  slowly  in  time  compared  to  the  symbol  time,  le., 
ak,(t)~akr(mT)  and  9k  r (r)  = 9k , (mT)  for  re  [mT,  (m  + l)T). 

If  we  want  to  achieve  spatial  diversity  for  fading  multipath  channels,  the  receiving 
antennas  must  be  spaced  sufficiently  far  apart  so  that  each  antenna  experiences  a different 
fading.  Usually  a physical  separation  of  at  least  10  wavelengths  is  required  between  two 
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antennas  in  order  that  each  antenna  obtains  a statistically  independent  fading  signal  [76]. 
In  this  case,  the  vector  channel  impulse  response  may  be  represented  as 


and  ° denotes  element-by-element  multiplication  which  is  called  Schur-Hadamard 
product  and  the  elements  of  the  Q.k  r(m)  are  complex  Lid.  (independently  identical- 

distributed)  Gaussian  random  processes  whose  amplitude  is  known  as  Rayleigh  fading. 

We  have  developed  a general  spatial  channel  model  that  includes  spatial 
information.  It  can  be  reduced  to  several  specific  cases  in  which  we  may  be  interested. 

In  a multipath  fading  channel,  if  the  multipath  delay  spread,  xm , defined  as  a 
maximum  propagation  difference  between  two  paths  is  larger  than  a small  fraction  of  one 
chip  interval,  the  channel  can  be  said  to  exhibit  frequency  selective  fading.  On  the  other 
hand,  if  the  delay  spread  is  much  smaller  than  one  chip  duration,  the  channel  can  be  said 
to  exhibit  frequency  nonselective  or  flat  fading. 

A frequency  selective  fading  channel  exploiting  spatial  diversity  with  multiple 
DO  As  per  user  where  rm  > Tc  is  used  as  our  general  channel  model.  The  general  channel 
model  of  (2.7)  can  now  be  reduced  to  five  special  cases  as  follows. 

• Flat  fading  channel  with  multiple  DO  As  per  user  (rm  « Tc ) 


(2.7) 


where 


(2.9) 


Additive  white  Gaussian  noise  (AWGN)  channel  with  multiple  DO  As  per  user 
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Mr,f)  = '£jei8t'<&(<pk,r)8(.T-Tk,)  (2.10) 

r=l 

Frequency  selective  fading  channel  with  a single  DOA  per  user  (zm  >TC) 


hk(z,t)  = <t>((pk)o 


£nv(<)5(f-r(,) 


r=l 


(2.11) 


Flat  fading  channel  with  a single  DOA  per  user  (rm  « Tc ) 
h*  (t,0  = (Q*  (f)  o <D  (<pk  ))5(  t-tk ) 
AWGN  channel  with  a single  DOA  per  user 

hk(r,t)  = e je>  <D  (<pk  )S( r-rk). 


(2.12) 


(2.13) 


The  received  signal  vector  at  a base  station  antenna  array  for  our  general  channel 
model,  r(r)  = [r(1)  (t)  r(2)(t ) r(P)  (r)]r  can  be  written  as 

r«)  = RejX  t te.,  (»)■ » )>.  <<  - }+"«).  C- 14) 


where  n(r)  = [rc(1)(0  nm(t ) n(P)(r)]r,  n(,,)(0  is  a white  Gaussian  noise 

waveform  with  two-sided  power  spectral  density  No/2  and  Re(X)  denotes  the  real  part  of 
X.  The  received  signal  vector  at  the  base  station  corresponding  to  the  channel  models  for 
five  special  cases  can  easily  be  obtained  from  the  general  received  vector  of  (2. 14). 


2.2  Discrete-Time  Complex  Vector  System  Model 


The  receiver  front-end  at  each  sensor  consists  of  a standard  IQ-mixing  stage 
followed  by  an  integrate-and-dump  section,  as  shown  in  Figure  2.4.  Ignoring  double 


frequency  terms,  the  equivalent  complex  received  sequence  at  each  sensor, 
r(p)  (/)  = r/p)  (/)  + jr^p)  (/)  can  be  expressed  as 
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K Kk 


r(,)  (/)  = Xi  C ( L(«  -1)  / W J H- J «-T„  )*  + »">(/) , (2.15) 

t=l  r=l  lc  ( ‘ 

where  with  |_xj,  we  denote  the  largest  integer  y such  that  y < x and 

p(kPr  (m)  = aff  (jnT) exp (jdff  imT)) exp(;2(p - 1 ){nd / A) sin (pkr for  ULA  and 

n(p)  (/)  is  a zero  mean  white  complex  Gaussian  sequence  with  variance 

<7 2 = £'[|/i(/’)  (/)|2|  = and  where  Eb  , = T is  the  transmitted  energy  per  bit  for 

1 J Tc  Eb.  1 


the  first  user. 

During  the  mth  bit  interval,  the  N chip  samples  are  accumulated  and  stored  in  a 
received  vector  at  the  pth  sensor, 

r<«(m)  = [r“()V),  r“(W-l), , if>(l)]e  C**.  (2.16) 


r^p\t) 


—►8) 

V2cos(wcr) 


► 

-V2sin(wcr) 


rTo* 

f:,(  v» 

t = ITC 


t = IT 


%\l) 


Figure  2.4  : Receiver  front-end  at  the  pth  sensor 
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The  equivalent  complex  received  sequence  vector  for  our  general  channel  model 
with  frequency  selective  fading  from  the  multiple  antenna  arrays  at  the  / th  chip  interval 
can  be  expressed  as 

r«  W = X £ (m>  ° WVx.r  Y (/ " T*.' )dt  + n* (/)  ’ (2- 17) 


t=l  r= 1 


where 


nm (/)  = [n®  (/),  ni2)(/), , < W • 


(*>  (i\tT 


(2.18) 


That  is,  the  / th  chip  samples  from  P antenna  sensors  can  be  stored  in  a received  vector, 
rM(Z)  = [/f(/),  C\l), , /fW  € . (2.19) 


The  equivalent  complex  received  sequence  at  the  base  station  corresponding  to  the 
channel  models  for  five  special  cases  can  easily  be  obtained  from  the  general  received 
vector  of  (2.17)  as  follows. 

• Equivalent  complex  received  sequence  with  multiple  DO  As  per  user  via  flat 
fading  (rm  «TC) 

r-  m = i % (a.,  (m)  ('-V  )*+”.(')  (2-20) 


*=1  r=l 


Equivalent  complex  received  sequence  with  multiple  DO  As  per  user  via  AWGN 


channel 


Equivalent  complex  received  sequence  with  a single  DOA  per  user  via  frequency 


selective  fading  (rm  > Tc ) 
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1 f (mN+l)T' 

sk 

T J(mN+l-l)T€  * 
1 c 


Y\ 

sk(t-Tk,)dt  +n  m(l)  (2.22) 


)) 


• Equivalent  complex  received  sequence  with  a single  DOA  per  user  via  flat  fading 


• Equivalent  complex  received  sequence  with  a single  DOA  per  user  via  an  AWGN 
channel 


For  the  fading  channel,  we  exploit  the  assumptions  that  the  fading  process  is 
independent  between  multiple  antenna  sensors  and  the  elements  of  the  Qt  r (m)  ° <t>(q>k  r ) 

in  (2.17)  are  complex  ii.d.  zero  mean  Gaussian  random  variables  with  unit  variance  if  the 
number  of  DO  As  of  the  incoming  signals  at  the  base  station  antenna  array  goes  to 
infinity.  In  order  to  show  that  the  elements  of  the  Qkr  (m)  ° 0(<pt  r ) are  complex  i.Ld. 

Gaussian-distributed  with  zero  mean  and  unit  variance,  we  consider  the  case  where  the 
number  of  DO  As  of  the  incoming  signals  at  the  base  station  is  large.  In  this  case,  we 
assume  that  sk  (t-tk  r ) ~ sk  (t  -tk ) . Thus,  (2. 17)  can  be  expressed  as 


«tc) 


s.it-T^dt  + nm  (/)  (2.23) 


(2.24) 


sk(t-zk)dt  + nm(l).  (2.25) 


Then,  define  Z(m)  = X (m)  + jY(m)  = £ (Q*,r  ("0  ° ®((Pk.r ))  where  X (m)  = Re(z(m)) 


and  Y (m)  = Im (z(m)),  where  Im(x)  denotes  the  imaginary  part  of  X.  If  the  number  of 
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paths  is  large  and  no  one  path  dominates  all  the  rest,  then  the  central  limit  theorem  can  be 


applied  to  both  X(m)  and  Y(m) , and  hence  both  can  be  well  approximated  as  Gaussian 


random  variables.  Assuming  that  the  second  moment  of  each  element  of  the  Z(m)  is 


equal  to  one 


E\z(>\m)f]  =1 


= 1 , the  elements  of  Z(m)  are  complex  Lid.  zero  mean 


Gaussian  random  variables  with  unit  variance  if  the  number  of  DO  As  of  the  incoming 
signals  at  the  base  station  array  goes  to  infinity.  In  this  case,  spatial  dependence  of  the 
impinging  signals  at  the  array  by  means  of  direction-dependent  phase  shift  between  any 
adjacent  antenna  elements  no  longer  exist  and  the  array  geometry  is  important  only  in 
that  it  provides  independent  fading  between  multiple  antenna  sensors.  That  is,  the  system 
model  for  fading  channels  contains  statistically  independent  spatial  samples  at  each 
sensor.  Then,  the  antenna  array  functions  as  a diversity  combiner  instead  of  as  a 
beamformer.  Therefore,  the  system  model  becomes  that  of  a diversity  combiner. 

That  is,  the  equivalent  complex  received  sequence  vector  in  (2.17)  for  diversity 
combining  may  be  represented  as 


Note  that  there  is  no  spatial  information  like  DOA  in  the  received  signal  for  the  diversity 
combining  whereas  DO  As  are  specified  for  the  adaptive  beamforming. 

Let  the  received  matrix  at  the  m th  bit  sampling  time,  r(m)  e , be  defined  as 


(2.26) 


r(m)  = [r„(A0  r.(W-l)  r„(l)] 


(2.27) 


or 


(2.28) 
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Let  the  contribution  of  the  k th  user’s  signal  to  r(p)(m)  be  r kp\m) . We  then  have 


r£p)(m)  = (a* (ru,Ttf2,-,TMj  )s[p\m))T , 


(2.29) 


where 


)—  [at.1^tj)  ai,2  fok,l  ) ai,l(T*.*t)  a*.2  2‘  (2.30) 

s^)(m)  = [4')(m)  s(kp}(m)  - (m)  s(p}Rt  (m)]T  eC2^.  (2.31) 

The  elements  of  skp)(m)  are  (»0  = fikP) (m)^2*-i (m)  and  = fikp?  (m)zlk(jn)  for 

r = l, 2,-  -,Rk,  where  z2*_i (m)  = (m) +dk(m- 1)]/ 2,  z2Jk (m)  = [rf* (m) - dk (m - 1)]/ 2 . 

The  columns  of  \k  are  the  contributions  from  the  k th  user’s  Rk  paths  to  r^im) . The 
form  of  a*! (zkr ) and  ak2{tkr)  are  defined  by  xkr  and  ck(l)  as 


au(T*.r)= 


5kr  ( 

-^P(/V+l,l)  + 


a*,2(T*,/-)- 


8k, r 


P(^r+1,-1)  + 


(2.32) 


(2.33) 


respectively,  where  rk,  = pkrTc  + 8kr  such  that  pkr  is  an  integer  and  8kr  e [o,  Tc)  and 
ck  e R**1  is  defined  as 


ck=[ck(N),ck(N-l), ,ct(l)]7 


The  permutation  matrix  P(r,a)e  R is  defined  in  block  form  as 


(2.34) 


P(r,a)  = 


0 IN.r 
alr  0 


a =±1, 


(2.35) 


where  I r denotes  the  rx  r identity  matrix. 
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The  received  vector,  r ip\m) , may  be  expressed  as 


r (p)  (m)  = (a s(p)  (m)  + n(p) (m))T  e ClxN , 


where 


A = [a,  A2  AjfJeR 


Nx 2R 


s(/’)(m)  = |^s1(p)r(/n)  s (P)T (m) 


sif>r(m)]reC2W 


and  where  R is  the  total  number  of  paths,  R = ^Rk . 


*= i 


By  assuming  that  the  first  user  is  the  desired  user,  the  received  vector, 
be  expressed  as 

r (p)  ( m ) = (Aj  s[p)  (m)  + e(p)  ( m)f , 

where  the  sum  of  the  MAI  and  the  additive  noise  is  denoted  by 

K 

e(p)  (m)  = Ats*p)  ( m ) + n(p)  (m) 


and  the  contribution  of  the  first  user’s  signal  to  r(p)(m)  can  be  written  as 

A 1s,(p)  (m)  = A;  (tu  ,t12  , • • • ,rlA  )b{p)  (m)z,  (m) , 

where 

)=[^u  ^1,2] 

^1,1  = [ai,i(ru ) au(Tu)  **'  3l,i (fi,*, )]  * i = l,2 


B1(p)(m)  = 


tiP\m)  0 
0 ti"\ m) 


tiP)  («)  = [tiP  (m)  tiP  (m)  • ••  /Jgj  (m)f , 


(2.36) 

(2.37) 

(2.38) 

(p\m),  can 

(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 


(2.45) 
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Zi(m)  = [zi(m)  z2(m)]re  {+1,0, -l}.  (2.46) 

2.2.1  Adaptive  Antenna  Array  or  Diversity  Combiner 

The  adaptive  antenna  array  or  diversity  combiner  system  consists  of  P identical 
sensors  with  assumptions  used  in  section  2.1  and  an  adaptive  signal  processor  that 
generates  a weight  vector  for  linearly  combining  the  array  outputs, 

wi  = [wi(1)  vvi(2) ]r  6CM,as  shown  in  Figure  2.5.  The  available  information  of 

the  desired  signal  is  used  to  calculate  the  tap  weights  of  the  adaptive  antenna  array  or 
diversity  combiner.  The  tap  weight  estimation  algorithms  of  the  adaptive  antenna  array 
will  be  covered  in  Chapter  3. 


t 


Available  information 


Figure  2.5  : Block  diagram  of  a P element  adaptive  antenna  array  or  diversity  combiner 
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2.3  Adaptive  Beamformer  or  Diversity  Combiner  Output  Data  Model 

The  adaptive  beamformer  output  during  the  mth  bit  interval  is  formed  by  a linear 
combination  of  the  received  signals  specified  by  the  tap  weights  of  the  adaptive  antenna 
array.  The  beamformer  output  of  the  received  signal  for  our  general  channel  model  with 
frequency  selective  fading  can  be  written  as 


where  each  element  of  y(m)  corresponding  to  the  l th  chip  interval  can  be  written  as 


y(ro)  = (wfr(»))T 


(2.47) 


=kw  ym(N-i)  y.d)]r6Cw, 


(2.48) 


where 


A*,r  ("t)  = nkr  (m)  ° &((pkr ) . 


(2.49) 


Now,  ym  (/)  can  also  be  expressed  as 


(2.50) 


where 


(2.51) 


and 


1 r (mN+l)T' 

J J (mN+l-l)Tt 


sk(t-rktr)dt+Yf"nm(l).  (2.52) 


Then, 


y(m)  = u(m)  + e(m) , 


(2.53) 
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where 

u(m)  = [um (N)  um(N-l)  um(l)J  (2.54) 

and 

e(m)  = [em (AO  e„(N  - 1)  em(l)J . (2.55) 

Thus,  u(m)  is  the  contribution  of  the  desired  user  and  can  be  written  as 

u(m)  = A xs[  (m) , (2.56) 

where 

s' (m)  = (m)  s[2(m)  •••  s[2R_y{m)  s[2Ki  (m)]r  e C2"1* . (2.57) 

The  elements  of  Sj (m)  are  s'l  2r_l(m)  = xltr(m)zi(m ) and  2r (m)  = xx >r (m)z2 (m)  for 
r = 1, 2,  •••,  R{ , where  xu(m ) = A1>r(m) . The  columns  of  A!  are  the 

contributions  from  the  first  user’s  Rx  paths  to  y(m) . For  instance,  if  dx (m)  = 
the  first  user’s  r th  path  will  contribute  the  vector  xu(m)d{ (m) au  (rl  r ).  On  the  other  hand, 
if  dx  (m)  = -dy  (m  - 1) , the  contribution  will  be  xUr  (m)dy  (m)a12  (rl  r ). 

Thus,  y(m)  can  be  written  as 

y(m)  = AjS,  (m)  + e(m)e  CNxl . (2.58) 

The  vector  e(m)  is  assumed  to  be  independent  of  the  desired  signal  component  at 
the  output  of  an  adaptive  antenna  array  and  to  be  circularly  symmetric  complex  Gaussian 
random  vector  with  zero-mean  and  arbitrary  covariance  matrix  Q that  satisfies 

£[e(m,)eH(/n;)]  =QSx(i-j),  (2.59) 

where  8K  is  the  Kronecker  delta.  The  unknown  covariance  matrix  Q models  both 
thermal  noise  and  all  other  interference  signals  including  MAI. 
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On  the  other  hand,  the  diversity  combiner  output  data  model  does  not  contain  the 
array  response  vector  and  hence  each  element  of  the  diversity  combiner  output,  y(m)  can 
be  written  as 


*(0  = X£ (m) -5- J sk ( t - Tkr )dt  + wf nM (/) . 


*=1  r= 1 


(2.60) 


CHAPTER  3 

TAP  WEIGHT  ESTIMATION  OF  THE  ADAPTIVE  ANTENNA  ARRAY  OR 

DIVERSITY  COMBINER 

This  chapter  investigates  beamforming  or  diversity  combining  techniques.  We 
consider  two  beamforming  techniques  to  calculate  tap  weights  of  an  adaptive  antenna 
array  for  ULA.  The  first  one  is  based  on  training  sequences.  Since  the  received  signals 
contain  the  propagation  time  delay  for  each  user’s  signal,  the  Direct  Matrix  Inversion 
(DMI)  algorithm  [42]  using  the  sample  correlation  matrix  and  sample  steering  vector  is 
modified  to  minimize  the  minimum  square  error  (MSE)  between  the  beamformer  or 
diversity  combiner  output  and  the  reference  signal  by  searching  N timing  phases  of  the 
reference  signal.  This  technique  is  also  used  for  diversity  combining.  We  will  look  at  the 
estimated  steering  vector  for  the  training  sequence  based  beamforming  under  simplifying 
assumptions.  The  other  one  is  a linearly  constrained  minimum  variance  (LCMV) 
beamformer  used  as  a performance  benchmark  of  the  first  one.  We  will  define  the 
performance  measure  for  the  adaptive  antenna  array  or  diversity  combiner  and  present 
numerical  results. 

3.1  Tap  Weight  Estimation  Algorithm 

We  consider  two  methods  to  estimate  the  tap  weights  of  an  adaptive  antenna  array. 
One  approach  is  a training  sequence  based  beamforming  using  the  MMSE  criterion.  The 
other  one  is  a DOA  based  beamforming.  Only  the  training  sequence  based  diversity 
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combining  is  used  for  diversity  combiner  because  the  diversity  combiner  does  not  process 
any  information  of  DO  A. 

3.1.1  Training  Sequence  based  Beamforming  or  Diversity  Combining  (Method  1) 

Given  the  known  training  sequence  dx  (m) , we  define  a reference  signal  for  the 
adaptive  beamformer  as 

t1(m,Q  = [fu(m,Q  r1>2(m,i)  •••  /uv(m,i')]=[c1  e RlxW,  (3.1) 

where  [x]*  means  components  a through  b of  x and  c,  °dj  is  the  element  by  element 

multiplication  of  Cj  and  dj  where  Cj  = [cf  cf c[  ] € Rlx("+1)Ar  and 

d,  = [dj (0)  dt(l) d^Mlje  R1*"*0"  with  d,(m)  = J,(m)«Ie  Rlx*  where  1 is  a 

row  vector  of  all  ones,  Le.,  d,  (m)  contains  the  N same  elements  obtained  by  N 
oversampling  the  mth  known  data  bit  and  the  phase  of  the  training  signal  could  be 
i = 0, 1, 2,  • • • , N - 1 . Then,  the  weight  vector  can  be  chosen  to  minimize  the  mean-square 

error  (MSE)  between  the  beamformer  output  wf  rm  (N  + l-n)  and  the  reference  signal 
tu  (m,  i) . Thus, 

eH (m,/)  = ru(m,i)- wfrm(A  + l-n),  n =1,2,  •••,#.  (3.2) 

Then, 

MSE(vf1,  i ) = \e„  (m,if  }=  1 - 2 Re{w?  }+  wfR^ , (3.3) 

where 

=£{rm(iV  + l-n)r"(A  + l-n)}  and  p„. = £:{rm(A  + l-n)/M(m,  /)}  . (3.4) 
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The  equation  for  the  weight  vector  that  minimizes  the  mean-square  error  is  well-known 
to  be 

ww  = • (3-5) 

In  order  to  get  i which  makes  the  MSE  to  be  smallest,  we  compute  the  MSE(i)  for  all 

1 M 1 M 

possible  i using  R„  = —7  Y, r(m)rH  (m)  and  prt , Y r(m)tf  (m,  i ) . Let  q be  the 

MN  “x  MN  m=1 

value  which  minimizes  the  MSE,  Le.,  q = arg  min;  MSE(i ) . With  the  Direct  Matrix 
Inversion  (DMI)  algorithm,  the  estimates  of  tap  weights  are  given  by 

wu=R;‘p^.  (3.6) 


3.1.2  Estimated  Steering  Vector  of  the  Training  Sequence  based  Beamformer 
Let  us  observe  the  estimated  steering  vector  for  the  training  sequence  based 
beamforming  (method  1)  over  an  AWGN  channel  Assuming  that  Rk  = 1 , rk  = 0 and 


<7  = 0, 


1 M 

Prt.9 

MN  m=i 


m= 1 (=1  [k=2  J 

^eJ6'em  + tLs  SCi  (N  ~ 1)rfi (m) Is  ■ ie> c* (m) + ”»2)  (/)l 


m=l  /=1 


(\)dx  (m)WjKeie'eKp-l)^  ck  ( l)dk  (m)  + n(mP)  (/)! 

M!S/  m= 1 1=1  [*=2  J 


(3.7) 


By  taking  the  statistical  average  of  p„  q , the  expected  value  of  the  estimated  steering 


vector  can  be  written  by 
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JpxeiB'e^ 

P n,q=  ^ . 

^PxeiB'ej(P~X)>h 
= JPxeiB'<b{(Pl). 

Therefore,  the  steering  vector  could  form  the  adaptive  antenna  array’s  beams  in  the 
direction  of  arrival  of  the  desired  signal. 


3.1.3  DOA  based  Beamforming  (Method  2) 

As  an  alternative  approach  to  the  adaptive  antenna  array  using  method  1,  we 
consider  a linearly  constrained  minimum  variance  (LCMV)  beamformer  [43]  which 
minimizes  the  output  power  subject  to  the  gain  in  the  direction  of  a desired  signal’s  r th 
path  being  constant,  Le., 

min^wfR^W!  subject  to  wf<I>(plr)  = 1,  (3.9) 


where  0(<pl  r ) is  called  array  transfer  vector  or  array  response  vector  of  the  r th  path  of 


the  first  user’s  signal  given  by 

<D(,plr)  = [i  e*1*  ej2*-'  

By  using  a Lagrange  Multiplier  method,  the  solution  to  (3.9)  is  given  by 

w RjfoOM 

Wl  <t>(<pu)uR  ~M<px,r)_ 


(3.10) 


(3.11) 


In  this  case,  side  information  which  is  the  DO  As  of  the  desired  user  r th  path,  <px  r , is 


required. 
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3.2  Performance  Measures  for  the  Adaptive  Antenna  Array  or  Diversity  Combiner 


To  evaluate  the  performance  of  the  adaptive  antenna  array  and  diversity  combiner, 
signal  to  interference  plus  noise  ratio  (SINR)  performance  is  considered.  The  input  SINR 
is  measured  at  the  first  element  of  the  antenna  array  whereas  the  output  SINR  is 
measured  using  the  weight  vector  estimated  by  the  training  sequence  based  beamforming 
(method  1)  and  DOA  based  beamforming  (method  2).  We  define  chip  sample  based 
signal  to  interference  plus  noise  ratios  as  following. 


where 


where 


Input  SINR  = ■ 

^(of  +|/i1)(/)|2 

(3.12) 

(3.13) 

MN  m=i  i=i 

(3.14) 

1 

MN  m=i  i=I 

(3.15) 

vi1)(/)  = 

— sAt-zlr)dt 

V " 1,r  T J (mN+l-l)Tt  1V  l'r 

(3.16) 

and 


1 r(mS+l)Tt 

T J (mN+l-l)T 

Ac 


Sk(t~Tk,)dt 


(3.17) 


and 


35 


Output  SINR  = 


fv.COl 


w 


jwf  n„(/)|  +|wffm(/)| 


(3.18) 


where 


v.W-k’W  v™(/)  vi"(()]r  (3.19) 

t.(l)  = [/■“(()  /“’«)  fi,n<r>Y-  (3.20) 


3.3  Tap  Wei£ht  Estimation  of  the  Adaptive  Antenna  Array  for  an  AWGN  Channel 

We  consider  the  beamformer  output  data  model  over  an  AWGN  channel  with  a 
single  DO  A per  user  of  (2.13)  which  is  one  of  our  five  special  cases  of  our  general 
model.  In  this  case,  y(m)  of  (2.58)  can  be  written  as 

y(m ) = D.z,  (m)  + e(m) 

(3 ' 

= -*1A1  (Tj  )z  , (m)  + e(m)  eC"*1, 
where  D,  =xl\1(r1)  and  jct  = *[l\eie'v/"Q>((px) . 


3.3.1  Numerical  Results 

The  simulated  system  uses  N=31  chips  per  bit.  The  code  sequences  are  Gold 
sequences.  The  performance  in  each  simulation  is  obtained  from  250  Monte-Carlo  runs. 
All  the  additive  noises  at  each  sensor  are  assumed  to  be  a zero-mean  white  Gaussian 
noise  with  power  spectral  density  of  No/2.  The  carrier  phases  9k , timing  offset  zk , DOA 
(pk  and  data  bits  of  all  users  are  independent  of  each  other.  Unless  explicitly  stated,  the 
same  conditions  and  assumptions  mentioned  above  will  be  used  for  our  MATLAB 
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simulations  throughout  this  dissertation.  The  near-far  ratio  (NFR)  is  defined  as 

Pk  / P{  for  k = 2,3, ,K . The  SINR  performance  of  an  adaptive  antenna  array  for  the 

AWGN  channel  in  (2.24)  with  Rk  = 1 is  evaluated  in  this  section  for  some  simulated 
scenarios.  We  use  the  ULA  where  the  distance  between  two  consecutive  sensors  for  the 
AWGN  channel  is  chosen  to  be  half  a carrier  wavelength  to  avoid  spatial  aliasing 
throughout  this  section. 

Figure  3.1  shows  the  average  input  and  output  SINR  with  randomly  generated 
DOAs  for  each  user’s  signal.  It  is  seen  that  the  average  output  SINR  is  increased  by  both 
methods  as  the  number  of  antenna  sensors  increases.  It  is  shown  that  method  1 and  2 
have  almost  the  same  average  SINR  performance  improvement  for  0 dB  and  10  dB  NFR 
and  that  method  1 achieves  less  increased  average  output  SINR  performance  as  compared 
to  that  of  method  2 for  NFR=20dB.  We  can  also  see  that  there  is  about  4 dB  SINR 
improvement  for  3 antennas  under  0 dB  and  10  dB  NFR. 

Figures  3.2,  3.3,  and  3.4  plot  the  SINR  with  randomly  generated  DOAs  for  each 
user’s  signal  versus  the  trial  number  for  the  different  NFRs,  respectively.  Method  1 and  2 
give  SINR  improvement,  but,  for  20  dB  NFR,  method  1 sometimes  causes  degradation  of 
SINR  due  to  poor  estimation  of  tap  weights  which  could  perturb  the  received  signal  even 
if  the  average  output  SINR  is  increased  for  20  dB  NFR  from  Figure  3.1.  So,  it  can  be 
expected  that  the  NFR  threshold  at  which  method  1 degrades  the  output  SINR  more  than 
one  trial  during  100  trials  exists  between  10  dB  NFR  and  20  dB  NFR. 

In  order  to  see  how  the  randomly  generated  DOA  of  each  user  affects  the  SINR 
performance,  especially,  for  the  20  dB  NFR  situation,  the  DOA  of  the  k th  user’s  signal  is 
assigned  by  the  expression,  DOA(k)  = 0 5x-kx  0.05* , which  implies  that  for  10  users 
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(K=10),  no  user’s  signal  has  the  same  DO  A of  any  other  and  hence  each  user  has  partial 
correlation  with  the  other  users  in  terms  of  DO  A.  From  Figure  3.5,  3.6  and  3.7,  we  can 
see  that  the  use  of  different  DOAs  for  each  user  gives  more  stable  output  SINR 
performance  improvement.  By  comparing  Figure  3.4  with  Figure  3.7,  it  is  obvious  that  if 
there  are  high  power  interfering  signals  from  undesired  users  with  the  almost  same  DOA 
as  that  of  a desired  user  due  to  the  random  generation  of  DOA  for  each  user,  the  adaptive 
antenna  array  can  not  suppress  them  and  more  degradation  of  SINR  occurs. 


Figure  3.1  : average  SINR  with  random  DOAs,  K=10,  Eb  l /N0= lOdB,  and  M=100 
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Figure  3.2  : SINR  with  random  DO  As,  NFR=0dB,  P=3,  K=10,  Eb , / iVo=10dB,  and 

M=100 
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Figure  3.3  : SINR  with  random  DOAs,  NFR=10dB,  P=3,  K=10,  Eb  x /JVo=10dB,  and 

M=100 
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Figure  3.4  : SINR  with  random  DOAs,  NFR=20dB,  P=3,  K=10,  Eb  l / Na= lOdB,  and 

M=100 
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Figure  3.5  : SINR  with  fixed  DOAs,  NFR=0dB,  P=3,  K=10,  Eb  l / No=10dB,  and  M=100 
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Figure  3.6  : SINR  with  fixed  DO  As,  NFR=10dB,  P=3,  K=10,  Ebl  /Ao=10dB,  and 

M=100 
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Figure  3.7  : SINR  with  fixed  DO  As,  NFR=20dB,  P=3,  K=10,  Eb , / Ao=10dB,  and 

M=100 
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3.3.2  Conclusions 

The  adaptive  antenna  array  can  suppress  the  MAI  by  adaptively  controlling  the 
antenna  beam  pattern  using  the  known  DOA  of  the  desired  signal  or  the  known  training 
sequence  even  if  the  desired  DS-CDMA  signal’s  DOA  is  unknown  and  hence  the 
adaptive  beamformer  output  has  less  interference.  We  have  seen  that  the  adaptive  antenna 
array  increases  the  average  SINR  performance  and  the  training  sequence  based 
beamforming  and  DOA  based  beamforming  achieve  almost  the  same  improved  SINR 
performance  for  0 dB  and  10  dB  NFR.  But,  when  method  1 is  used,  there  are  some 
degradations  of  SINR  performance  for  20  dB  NFR.  We  will  see  how  these  observations 
of  SINR  performance  are  related  to  the  performance  of  the  propagation  delay  estimation 
problem  in  chapter  5. 

3.4  Tap  Weight  Estimation  of  the  Diversity  Combiner  for  Fading  Channels 
3.4.1  Numerical  Results 

The  SINR  performance  of  a diversity  combiner  for  a time- varying  fading  channel  is 
evaluated  by  assuming  a single  path  ( Rk  = l)  flat  fading.  In  the  simulations  throughout 
this  dissertation,  the  fading  process  is  generated  by  passing  white  gaussian  noise  through 
a third  order  low  pass  filter.  In  this  simulation,  the  data  rate  is  assumed  to  be  9600  bps 
and  Eb  l / N0  =10dB  is  used. 

As  the  number  of  antenna  sensors  increases,  the  average  output  SINR  performance 
is  shown  to  increase  in  Figure  3.9  when  we  use  M=50  (observation  bits),  K=20  (users), 
NFR=0dB  and  fd  (Doppler  frequency)  = 70  Hz  (or  a normalized  doppler  rate  of  0.00729) 
where  a Doppler  rate  of  70  Hz  would  result  if  the  system  had  a carrier  frequency  of  900 
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MHz  and  a transmitter-receiver  relative  speed  of  84  km/h.  We  can  see  in  Figure  3.10  that 
when  M=50,  P=4,  K=20  and  NFR=0dB  are  used,  the  amount  of  the  output  SINR 
improvement  decreases  as  the  Doppler  frequency  increases. 

3.4.2  Conclusions 

The  diversity  combiner  can  suppress  the  interfering  signals  by  exploiting  the  spatial 
diversity  of  the  received  signals.  It  is  seen  that  the  diversity  combiner  can  improve  the 
average  output  SINR  performance  for  0 dB  NFR  and  the  output  SINR  improvement  is 
reduced  as  the  fading  rate  increases. 
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Figure  3.10  : average  SINR  with  Eb  x / Afo=10dB,  NFR=0dB,  M=50,  P=4  and  K=20 


CHAPTER  4 

CRAMER-RAO  BOUND  (CRB)  FOR  PARAMETER  ESTIMATION 

The  main  contribution  of  this  chapter  is  the  derivation  of  the  CRB  for  parameter 
estimation  in  a DS-CDMA  system  on  an  AWGN  channel  and  fading  channels  when  we 
use  multiple  antenna  array  sensors.  This  CRB  is  used  as  a performance  bound  for  the 
accuracy  of  the  propagation  delay  estimates. 

4.1  CRB  for  an  AWGN  Model 

We  will  consider  two  kinds  of  the  Cramer-Rao  bounds  for  the  AWGN  channel.  The 
first  one  is  the  CRB  for  the  parameter  estimates  of  the  beamformer  output  data  model 
after  linearly  combining  the  received  data  from  P antenna  sensors.  This  CRB  is  a function 
of  tap  weights  of  the  adaptive  beamformer  and  is  derived  from  the  beamformer  output 
data  even  if  the  original  data  are  received  from  P antenna  elements.  That  is,  this  CRB  is 
the  lowest  bound  on  error  covariance  of  the  estimator  of  a parameter  when  we  use  the 
beamformer  output  data  model  in  (3.21)  and  is  referred  to  as  post-combining  CRB.  Since 
this  CRB  is  affected  by  the  adaptive  beamforming  process,  we  consider  another  bound  on 
error  covariance  of  any  unbiased  estimator  of  a parameter  with  the  received  data  model 
from  P antenna  sensors.  So,  the  second  one  is  the  CRB  regardless  of  the  choice  of 
combining  process  which  is  referred  to  as  pre-combining  CRB.  That  is,  this  CRB  is  the 
best  possible  performance  for  any  unbiased  estimator.  By  comparing  the  first  CRB  (post- 
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combining  CRB)  with  the  second  one  (pre-combining  CRB),  we  examine  how  much 
information  about  the  signal  timing  is  lost  in  the  process  of  beamforming. 

4.1.1  Post-combining  CRB 

Since  we  are  applying  the  LSML  timing  estimator  to  the  output  data  of  an  adaptive 
antenna  array  after  the  received  data  are  combined,  we  describe  the  CRB  for  the 
parameter  estimates  of  the  data  model  for  the  output  of  an  adaptive  antenna  array 
presented  in  (3.21). 

If  we  consider  the  parameter  as  the  propagation  delay  r, , the  carrier  amplitude 
7i  and  phase  dx  of  the  desired  user’s  signal,  the  parameter  vector  can  be  written  as 

r\  = t7i  ^]r  • The  extended  Bangs’  formula  for  the  ijth  element  of  the  Fisher 
information  matrix  (FIM),  J,  is  of  the  form  [44,48] 

{J  % =/r(Q-1V1QQ-1VJ.Q)+2Re[V1.(zfDa)Q-1V>(l>z1)]  , (4.1) 

where  V,X  denotes  the  gradient  of  X with  respect  to  the  ith  unknown  of  the  likelihood 
function,  tr(X)  denotes  the  trace  of  X.  Note  that  the  FIM,  J is  block  diagonal  since  Q 
does  not  depend  on  the  parameters  in  (Dz! ),  and  (Dz, ) does  not  depend  on  the  elements 
of  Q. 

The  submatrix  of  the  Fisher  information  matrix  corresponding  to  the  parameter 
space  is 


1 

0 

J= 

0 

J6A 

0 

0 

A*. 

where 
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J7l7l  (m)AfQ  ‘AjijCw)}  (4.3) 

m=l 

Jm  =2X|x,|{.f(m)A"Q-,A^,(m)}  (4.4) 

m= 1 

JsA  =2X|x1|2tf,(/n)A1ffQ-1A1zI(m)}  (4.5) 

m- 1 

=2^|a,|2^|H  (m)A^Q'lAfz1(m)}  (4.6) 

m— 1 

au  a,  2 eR^with 

au,/  = ^-— = ^[p(Pi  +l,l)-P(pl,l)]c1  (4.7) 

oTi  Tc 

and 

+ l,-l)-F(f>„-l)]c,.  (4.8) 

&l\  Tc 

Then,  the  CRB  for  r\  is  given  by 

CRB(77)  = J"1.  (4.9) 

Note  that  the  post-combining  CRB  can  be  derived  in  the  exact  same  way  as  in  [17]  and  it 
is  the  same  as  the  CRB  in  [17]  except  for  x, . Remember  that 

A,=V^'«"'wf<I>(«>1)-  (4.10) 


4.1.2  Pre-comhining  CRB 

We  derive  the  Cramer-Rao  bound  for  the  AWGN  channel  that  M received  vectors 
from  P antenna  sensors  can  get  no  matter  what  estimators  or  beamformers  are  used.  In 
order  to  derive  the  CRB  for  the  AWGN  channel  given  M received  vectors  from  P antenna 
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sensors  without  considering  tap  weights  of  an  adaptive  antenna  array,  we  follow  the 
procedure  in  [45]  and  define  the  received  vector  in  a little  different  manner  from  the 
previous  one.  We  assume  that  Rk  = 1 . 


Let  the  received  vector  during  the  m th  bit  interval,  r(m)  e CNM 
vector  n(m)e  CNPxl  be  defined  as 

r(m)  = [r(1)(m)  r(2)(m)  •••  r(/>)(m)  ]r 

and 

n(m)  = [n(1)(m)  n(2)(m)  •••  n(/>)(m)  ]r, 

where 

""’(m)  = k"W  - ni'»(l)]  € C1** 

The  r(m)  can  be  written  in  a more  compact  form  as 

K 

f (m)  = fk  ( m ) + n(m)  = B(t)E  (y,  0)z(m)  + n(m) , 

k= 1 

where  r*  (m)  is  the  contribution  from  the  k th  user  to  r(m)  and 

B(r)  = [b,  b2  - b2K]eCNP*2K, 

E = C2™\ 

z(m)  = [z,  (m)  z2  (m)  • • • z2K  (m)f  e {+ 1,  0,  - 1}2™ , 

where 

t = [ti  t2  xj 

r = Iri  r2  •••  YkJ 

o = [el  e2  ...  eKJ 


and  the  noise 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


(4.20) 
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where  = ■s[l\ej8> , yk  = -JK  and 

•(*»)  T JirccLj-y-l  in <x»fyt-(f 

a*.,(Tt)r  e ^ W7*)  •••  <?  1 Ja*,i(T*> 


h2t-i  ~ 


]KQX*\{ 


6 C 


WPxl 


(4.21) 


and 


h2*  - 


J* «(»,)  ) j«co»fyt  (P  ‘)2g  1 _ 

a*,2(T*)r  e ^ '«*.»(** ) •••«''  Jat>2(rt) 


eC 


A7>xl 


(4.22) 

The  parameter  estimation  problem  can  be  formulated  as  follows:  Given  M 
observations  of  the  received  vector  r(m)  from  an  antenna  array,  estimate  the  parameter 

vector  ¥ = fr2  yT  dT  rT  ]r  e R3Ar+1  , 

The  CRB  is  a lower  bound  on  the  error  covariance  matrix  of  any  unbiased  estimator 

A 

of  a parameter  . The  covariance  matrix  of  the  unbiased  estimator  *¥  satisfies 

= CRB('T) , (4.23) 


JW- 

f 

E 

f d In  L(r)  Y d In 

t‘ 

y 

< 

.1  9t  A **  J 

> 

where  Je  R(1+3a:)x(1+3a:)  is  Fisher  information  matrix  (FIM). 

Since  the  derivation  of  the  CRB  is  somewhat  lengthy,  we  defer  the  details  to 
Appendix  A.  1.  It  is  shown  in  Appendix  A.1  that  the  Fisher  Information  matrix  (FIM), 
Je  r (1+3X)x(i+3 jc)  may  written  as 

' MNP 


J = 


a 

0 

0 

0 


0 0 0 

J?7  ^y0  Jjt 

J ye  Jee  Jfo 

Tr  Tr  T 

J)T  J9r  JTT 


(4.24) 
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where  the  matrices  J,,,,  J^,  J^,  Jw,  J<*,  e RM  are  defined  as 


J„  =^-XRe{ztf(m)E«BHBErZ(m)} 

G m-l 

(4.25) 

G m= 1 

(4.26) 

= ^jrRe{z"(ro)E?B*GEZ(m)} 

O'  m=\ 

(4.27) 

J*  = ^TfJR^H(m)E^BHBEgZ(m)} 

G m=i 

(4.28) 

J st  = ^J,R^\m)EHgBHGEZ{m)} 

G m=l 

(4.29) 

J xt  = ^YfJRe{zH(m)EHGHGEZ(m)} 

(4.30) 

G “ 


where 


Er  = diag(exp(jdl ), exp(jOl ), ,exp(jGK),exp(jdK))e  C2Kx2K , (4.31) 

Ee  =yEeC2Jfx2Jf  (4.32) 

and  where 

G = [g,  S,  - blsC*™  (4.33) 

is  defined  as 


_d*>2k-l 

dzk  ’ 


82*-!  — ll  » 82*  - 


_ d b2k 


dxk 


(4.34) 


and  where  Z(m)  g {- 1, 0,  l}2*x*  is  defined  as 


(Z  (m)\i,k)  = 


Zu-i  ("0  i/  1 = 2*  - 1 
Z2i  (w)  if  i = 2k 
0 otherwise 


(4.35) 
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Note  that  when  P = 1 , the  pre-combining  CRB  is  the  same  as  that  obtained  in  [20]  and 
that  when  Q = tr2I , the  pre-combining  CRB  is  the  same  as  the  post-combining  CRB  if 
there  is  no  tap  weight  for  P = 1 and  K = 1 . 


It  can  be  shown  by  following  similar  calculations  used  in  [20]  that  the  asymptotic 
Fisher  information  matrix,  J„,  is 


J„  = 


lim  J = 

M — *<*> 


MNP 

a4 

0 

0 

0 


0 

0 
T T 


0 

0 

Jee,. 

0 


0 

?*.■ 

0 


(4.36) 


where  Jn „ eRM  is 


•>».-  = Um  A|;Re{Z»(m)E»B''BErZ(m)} 


2 

O7  m=l 

= Re{Lw  diag  (b  H B 


(4.37) 


since  E^E r = I2K  and  LH 1 2KL  = 21 K.  Furthermore,  since  E"E  = |E|,  EfEe  =|e|2  and 
EhE  = |E|2  where  we  by  |E|  denote  the  matrix  whose  elements  are  |e|(i,  j)  =|E(i,  j) | , 


J w. o » J ee.oo i J tt.oo  ^ R are 


i jt,~  > J ee,«  > « rr,1 


Jn-  = 1 = ^-Re{L"d;as(B*G^iL"lElL} 

= Re{Lw^  (bhG^Vw  } 
J».-  = Km  Jw  =^-Re{Lff^(B//B>Lw} 

Af  — ►°° 

JCT„  = lim  = ~Re{LH dag  (GHG)LW} 

A/  >°°  ^ 


(4.38) 


(4.39) 


(4.40) 
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where  Le  {0,  l}2ArxA:  and  the  power  matrix  WeRw  are  defined  as 


L(i,*)  = 


1 

1 

0 


if  i = 2k-l 

if  i = 2k 

(4.41) 

otherwise 

.f,}. 

(4.42) 

We  can  invert  the  asymptotic  Fisher  information  matrix  analytically  since  all  blocks  are 
diagonal.  Moreover,  the  structure  of  J„  is  such  that  the  bound  on  the  variance  of  an 
unbiased  estimator  of  yx , 6X,  or  r,  is  independent  of  the  parameters  for  any  other  user. 
Especially,  the  bound  on  the  variance  of  a r,  estimator  is  independent  of  Pk  for 
it  = 2, 3,  • • ',  K . Since  our  interest  is  in  the  parameter  vector  r , the  CRB  of  r by  using 
(A.29)  can  be  obtained  as  follows: 

0 


CRB_1(t)  = 


CRB~\tx) 

0 

0 


CRB-\z2) 


0 


0 

0 

CRB'1  (tk  )j 


€ R 


KxK 


(4.43) 


where 


CRB  (rk  ) - 2 {(g2*-i§2*-i  "*"82*82* )+  [^e (b 2*— 1 8 2*-i  + b^g2* )]  (b^t-ib^-i  + b2kb2k )"  }. 


(4.44) 

Therefore,  the  CRB  using  the  asymptotic  Fisher  information  matrix  is  near-far  resistant 
and  independent  of  the  number  of  users  on  an  AWGN  channel.  Furthermore,  as  indicated 
in  subsection  4.1.3,  the  actual  CRB  is  very  close  to  the  asymptotic  expression  for  rather 


modest  values  of  MP . 
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4.1.3  Numerical  Results 

The  performance  in  each  simulation  is  obtained  from  10,000  Monte-Carlo  runs. 
Figure  4.1  shows  the  comparison  between  the  pre-combining  CRB  and  post-combining 
CRB  when  methods  1 and  2 are  used  as  a function  of  the  number  of  observation  bits.  It  is 
shown  that  both  the  pre-combining  CRB  and  the  post-combining  CRB  decrease  as  the 
number  of  observation  bits  increases  and  that  the  post-combining  CRB  depends  on  the 
NFR.  It  is  seen  that  the  post-combining  CRB  when  the  method  2 is  used  in  the  0 dB  NFR 
situation  does  not  have  much  difference  from  the  pre-combining  CRB  from  Figure  4.1  as 
well  as  Figure  5.2  and  5.3  in  chapter  5.  The  fact  that  the  post-combining  CRB  is 
approaching  the  pre-combining  CRB  tells  us  that  the  beamformer  does  not  destroy  much 
information  about  the  timing  of  the  signal  and  hence  the  linear  combination  of  the 
received  signal  from  P antenna  sensors  using  the  tap  weights  with  method  1 or  2 and  then 
using  the  LSML  estimator  could  be  a good  propagation  delay  estimation  method  under 
the  condition  of  the  reasonable  power  control  Although  the  asymptotic  FIM  of  (4.36)  can 
be  used  to  calculate  the  performance  bound  because  the  finite  sample  CRB  approaches 
the  asymptotic  CRB  for  MP  > 20 , our  simulations  were  performed  using  the  finite 
sample  CRB  of  (4.24). 

4.1.4  Conclusions 

The  fact  that  the  post-combining  CRB  determined  by  two  methods  used  in  Chapter 
3 for  calculating  the  tap  weight  vector  of  the  adaptive  beamformer  is  not  much  different 
from  pre-combining  CRB  under  the  perfect  power  control  condition  tells  us  that  the 
adaptive  antenna  array  plus  LSML  timing  estimator  could  be  a good  timing  estimator. 
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More  comparisons  between  the  post-combining  CRB  and  pre-combining  CRB  in  various 
scenarios  will  be  given  in  Chapter  5 when  dealing  with  the  timing  estimation  problem. 


Figure  4.1  : Post-combining  CRB  with  random  DOAs,  P=3,  K=10,  and  Eb , / Ao=10dB 


4.2  CRB  for  Fading  Channels 


In  this  section,  the  pre-combining  CRB  from  the  subsection  4.1.2  is  generalized  to 
fading  channels.  Our  approach  to  the  extension  to  fading  channels  is  to  consider  members 
of  the  set 


{$?("*)  : m = \,2,-,M-k  = \,2,-,K-p  = \,2,-,P-r  = \,2,-,Rk  } (4.45) 
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to  be  unknown  and  deterministic  parameters.  The  parameter  vector  is  then  formed  as  the 
z vector  augmented  with  <7 2 and  the  elements  of  fe?  ("»)}• 

We  derive  the  CRB  for  fading  channels  given  M received  vectors  from  P antenna 
sensors  by  following  the  similar  procedure  in  [45].  In  order  to  derive  the  CRB  for  fading 
channels,  we  define  the  received  vector,  r(m)  as  in  (4.11)  where  the  equivalent  complex 
received  sequence  vector  for  diversity  combining  of  (2.26)  is  used,  Le.,  each  element  of 

a»,(m)Vn  is  = • 

After  some  straightforward  calculations,  we  can  formulate  the  contribution  from  the 
kth  user  to  r(m)  as  rt  (m)  = XkEk  (m) , where 


Xk  diag{Ak,Ak,~ • 

Aje 

(4.46) 

St(ro)  = [s!t1)r(m)  E[2)T(m) 



(4.47) 

where  Ak  is  defined  as  in  (2.30)  and 

( m ) = \fikx  ("QZu- 1 (»»), 1 3ft ( m)z2k ( m ), 

— . P(kPl  ('«)z2i_1  (m),  ( m)z2k (m)]F  eC2^. 

(4.48) 

The  received  vector,  r(m) , may  be  expressed  as 

r(m)  = XS(m)  + n(m) , 

(4.49) 

where 

x = [x,  x2  

XK]e  RNPx2RP 

(4.50) 

S(m)  = [H1r(m)  E2T  (m)  • 

EJcr(m)]r  e C1*™ 

(4.51) 

and  where  R is  total  number  of  paths 
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K 


«=£«*• 


(4.52) 


4.2.1  CRB 

The  derivation  of  the  Cramer-Rao  bound  (CRB)  is  somewhat  lengthy,  and  we 
therefore  defer  the  details  to  the  Appendix  A.2. 

The  basic  assumption  made  for  the  AWGN  case  in  subsection  4.1.2  are  made  here 
as  well : a2  and  zk  r are  assumed  to  be  unknown  and  deterministic  parameters  that  are 
fixed  during  the  observation  time  interval  We  assume  that  rk/.  # pTc  for  an  integer  p, 
which  guarantees  that  X is  differentiable  with  respect  to  zkr . Furthermore,  we  assume 
that  \fikp,  (m)}  is  unknown  and  deterministic  for  m = 1, 2,  •••,  M , k = 1, 2,  •••,  K , 
p = l,  2, -,P  and  r = l,2,  •••,/?*. 

The  real  part  of  a complex  quantity  x is  denoted  by  x and  the  imaginary  part  is 
denoted  by  x . The  parameter  vector  4*  e R1+2™+*  is  defined  as 


= PT  tt}\ 


(4.53) 


where 


IRPMxl 


(4.54) 


(4.55) 


(4.56) 


(4.57) 


and 
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-k 


T 


<r e r 


«xl 


— \pk,l  ^k,2  f € R ' * 

The  Fisher  information  matrix  (FIM),  J e is  found  as 


(4.58) 

(4.59) 


J = E 


^ lnL(r)  Y <9  lnL(f)Y 

, & A & > 


(4.60) 


where  In  L(r)  is  the  log-likelihood  function  (conditioned  on 


z = ^zr(l)  zr(2)  •••  i (M)j  where  z(m)  = [z1(m)  z2(m)  •••  z2K(m)Y  e 

{+ 1, 0,  - 1}2™ ) of  r = [ fT  (1)  r r (2)  • • • fT  (M )]  with  respect  to  *F . Thus,  the 
resulting  CRB  should  be  interpreted  as  being  conditioned  on  both  the  data  sequence  and 
the  fading  processes. 

It  is  easy  to  show  that  the  likelihood  function  conditioned  on  z is 


( 1 M 


"TlXll  r(m)-XH(m)||2 

° m=l 


(4.61) 


and  ignoring  constants  independent  of  , the  log-likelihood  function  is 

1 m 

lnL(r)  = -MA/Plna2 r£||  r(m)-XS(m)||  . 


O m= 1 


(4.63) 


In  Appendix  A.2,  it  is  shown  that 


<91nL(r)  MNP  , 1 ^ 


m= 1 


^-lARe{z«(m)X»a(m)} 

MU4.wx-iwj 


(4.64) 


(4.65) 


(4.66) 
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dlnL(r)=^|.Re^„  (m)H»fi(m)} 
dT  O1 

(4.67) 

where  H€  rNFx2RP  is  defined  as 

H = [H,  H2  

H,] 

(4.68) 

H*  =diag(\Tk,ATt',— 

(4.69) 

^ _ 3A* 

(4.70) 

r*' 

and  where  Zk(m)e  {-1, 0,  l}2S‘  *R" , Z(m)  e {- 1, 0,  \yRP*RP , 

rt(m)e  C2/w  and 

r(m)e  C2*™  are  defined  as 


Z(m)  = diag  (Dt  (m),  D2  (m), 

•»Dj f(m)) 

(4.71) 

•3 

n 

✓ — \ 

6 

'W 

•w 

Q 

(m),Zk(m),-- 

Zt(m)) 

(4.72) 

*24-1 0»)  1/ 

i = 2 j-l 

(Z 

= - 

z2*(™) 

i = 2j 

(4.73) 

0 otherwise 

r(m)  = diag(rx(m),  T2(m), •••• 

S,(m)) 

(4.74) 

Phpy(m)z2k-i(m) 

if  i = 2 j-l 

(r*  j)  = • 

P$(m)z2k(m) 

if  i = 2; 

(4.75) 

0 

otherwise 

where  p - 


i-l 

2Rk_ 


+ 1 and  r = j-nxRk  where  n is  the  integer  part  of  the  quotient. 


j~  1 
Rk  ‘ 

The  FIM,  J,  can  now  be  formed  by  taking  the  expected  value  of  all  possible  outer 


products  of  the  gradients  (4.64)-(4.67).  After  some  manipulation  we  conclude  that 
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J = 


MNP 

cr4 

0 

0 


0 


\T 

'fix 


0 

hr 

Jrr 


(4.76) 


where 


3 pp  — diagij  pp  (1),  J pp  (2), , 3 pp  ( M ))  (4.77) 


J (m)  = — 


Re{zw(m)XffXZ(m)} 

- Imr  {zH  (m)XHXZ(m)} 


- Im{zff  (m)XHXZ(m)} 
Re{z"(m)X*XZ(m)} 


(4.78) 


and 


3pr=[3Tpt(l)  3 Tpt (2)  3Tpx{M)]T 

Re{z"(m)X*HT(m)} 
Im{ztf(m)X"Hr(m)l 

Jrr=4-XRe[rW(m)HflHr(m)]. 


(4.79) 

(4.80) 

(4.81) 


In  order  to  get  the  CRB  for  an  unbiased  estimator  of  z , we  invert  the  FIM 


J‘  = 


— — 0 
MNP 

0 x 


0 

x 


0 x CRB(z) 


(4.82) 


where  CRB(z)e  R***  i s the  lower  bound  right-hand  RxR  block  matrix  of  J ‘.In 
(4.82),  x denotes  matrices  that  are  of  no  particular  interest  to  us.  It  is  shown  in  Section 
A.  2 that 

, (4.83) 

where  f is  an  unbiased  estimator  of  r and  then 


CRB~l  (r)  = 4-ERe[r"  (m)HwP^(m)HT(m)] 
<7  m=l 

where  PxZ(m)  is  the  projection  matrix  onto  [range  (XZ(m))]J' , 


Px1,,.,  = /»,  - XZ(m)(z*  (m)X"XZ(m)r  Z"  (m)X"  . 

Numerical  results  of  the  CRB  for  fading  channels  will  be  given  in  Chapter  5 when 
dealing  with  the  timing  estimation  problem. 
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(4.84) 


(4.85) 


CHAPTER  5 
TIMING  ESTIMATION 


In  this  chapter,  the  LSML  timing  estimator  which  follows  the  adaptive  antenna 
array  and  diversity  combiner  is  introduced  for  estimating  the  propagation  time  delay  of  a 
desired  signal.  The  CRBs  derived  in  Chapter  4 are  used  as  a performance  lower  bound. 
By  comparing  two  kinds  of  CRB,  we  will  investigate  the  performance  of  our  approach, 
Le.,  the  adaptive  beamformer  or  diversity  combiner  plus  LSML  timing  estimator. 

Figure  5.1  is  a block  diagram  of  an  adaptive  antenna  array  or  diversity  combiner 
plus  LSML  timing  estimator  for  estimating  the  propagation  time  delay  of  the  desired 
user’s  signal.  The  received  signal  from  P antenna  sensors  is  stored  over  the  observation 
time  interval  [0,  MT),  and  then  the  tap  weights  of  the  adaptive  antenna  array  or  diversity 
combiner  are  computed.  The  beamformer  output  which  is  formed  by  linearly  combining 
the  array  output  has  less  interference.  Then,  the  LSML  timing  estimator  is  applied  to  the 
beamformer  or  diversity  combiner  output  in  order  to  get  the  propagation  delay  estimate 
of  the  desired  user’s  signal. 

We  present  the  MATLAB  simulation  results  for  various  scenarios.  We  will  evaluate 
the  performance  of  our  approach  on  an  AWGN  channel  and  then  on  a fading  channel, 
respectively.  For  the  AWGN  channel,  the  performance  of  the  adaptive  beamformer  based 
LSML  timing  estimator  is  compared  with  that  of  a single  antenna  based  LSML  timing 
estimator.  Even  if  the  MASE  algorithm  in  [26]  provides  good  performance  on  the  flat 
fading  channel  with  receiver  diversity,  it  is  unknown  whether  the  MASE  algorithm  still 
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works  well  on  a channel  including  direction  of  arrival  (DO A)  information  and/or  a 
system  without  receiver  diversity.  Therefore,  the  MASE  algorithm  is  applied  to  an 
AWGN  channel  where  the  only  different  information  between  any  antenna  sensors  is  a 
DOA-dependent  phase-shift  and  its  performance  is  compared  with  that  of  the  adaptive 
antenna  array  plus  LSML  timing  estimator.  On  the  other  hand,  the  LSML  timing 
estimator  based  on  the  diversity  combiner  which  can  reduce  the  effects  of  channel  fading 
on  a desired  signal  is  applied  to  flat  fading  channel  with  a single  path  and  the 
performance  of  the  MASE  algorithm  developed  for  a single  path  flat  fading  channel  is 
compared  with  that  of  our  approach.  The  MASE  algorithm  considers  only  the  case  of  a 
single  path  flat  fading  channel  and  do  not  consider  the  case  of  frequency  selective  fading 
due  to  multipath  propagation.  Therefore,  we  extend  the  MASE  algorithm  to  a frequency 
selective  fading  channel  for  joint  estimation  of  multipath  timings  of  a desired  user’s 
signal.  And,  the  MUSIC  timing  estimator  in  [47]  is  also  extended  to  the  case  of  using  a 
multiple  antenna  array  for  estimating  multipath  timings  in  frequency  selective  fading 
channels. 


5.1  Timing  Estimation  for  an  AWGN  Channel 

The  following  subsection  deals  with  the  propagation  delay  estimation  problem 
based  on  an  adaptive  beamformer  in  a DS-CDMA  system  over  an  AWGN  channel.  The 
beamformer  output  data  model  in  (3.21)  is  used. 
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Adaptive  Beamformer  or  Diversity  Combiner 


Figure  5.1  : Block  diagram  of  an  adaptive  antenna  array  or  diversity  combiner  plus 

LSML  timing  estimator 


5.1.1  LSML  Timing  Estimator 

The  LSML  timing  estimator  originally  developed  on  the  basis  of  a single  antenna 
[17]  can  be  derived  from  the  adaptive  beamformer  output  data  model  of  (3.21)  and  hence 
is  given  by 


A 


= arg  max{Ti) 


||®  u (Ti  + ®u  (Ti  )^2 1| 

||®u(Ti)||2  +||®i.2(ti)||2 


(5.1) 


where 


di  =Q_1/2dj,  d2  = Q-1/2d 


(5.2) 
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and 


5u(ti)  = Q 1/2aM  (Ti),  au(T,)  = Q l/zau(T t ) 


-l/2„ 


(5.3) 


where  dj  and  d2  denote  the  first  and  the  second  columns  of  D,  respectively,  with 


and 


where 


D^R;1  (5.4) 


Q=K-K*tK 


(5.5) 


M m=l 


(5.6) 


A.-S-t-.ttK'W 

^ m=l 


(5.7) 


M m=i 


(5.8) 


5.1.2  MASE  Algorithm 

The  MASE  algorithm  [26]  is  an  asymptotic  (for  large  number  of  data  samples)  ML 
estimator  for  a multiple  antenna  sensor  array  which  exploits  receiver  diversity  on  the  flat 
fading  channel  The  MASE  algorithm  is  given  by 


Tj  =argmax{Ti)- 


[aufrild^  +a1g2(r1)d1(/,)| 

iMor+MMf 


(5.9) 


where 


d{p)  = Q_1/2d[p),  d\p)  = Q-1/2d^ 

*u  (*i ) = Q"1/2au (Ti  )>  »1,2  (Tj ) = Q-1/2a12  (t,  ) 


(5.10) 

(5.11) 
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D(,)  =R^R“‘  =[dif)  dH 

(5.12) 

r />=! 

(5.13) 

Ry, 

M m=l 

(5.14) 

R ,,,,  4ir(f)r(m)rW*(m). 

M m=l 

(5.15) 

5.1.3  Numerical  Results 

The  performance  of  an  adaptive  antenna  array  based  LSML  timing  estimator  over 
an  AWGN  channel  is  evaluated  in  this  subsection  for  some  simulated  scenarios.  We  also 
compare  the  adaptive  antenna  array  based  LSML  timing  estimator  with  the  single  antenna 
based  LSML  timing  estimator.  We  consider  a UCA  where  the  radius  r is  chosen  to  be 
half  a carrier  wavelength.  The  simulated  system  uses  N=31  chips  per  bit  and  Tc=  1 
operating  over  an  AWGN  channel.  All  the  additive  noises  at  each  sensor  are  assumed  to 
be  zero-mean  white  Gaussian  noise  with  power  spectral  density  of  No/2.  The  timing 
offset  Tk , DOA  (pk  and  data  bits  of  all  users  are  independent  of  each  other.  The 
performance  in  each  simulation  is  obtained  from  10,000  Monte-Carlo  runs.  In  our 
simulations,  the  same  known  training  sequences  were  used  in  calculating  the  tap  weights 
of  the  adaptive  antenna  array  with  DMI  algorithm  and  in  estimating  the  propagation 
delay  of  the  desired  user  in  the  LSML  estimator. 

We  consider  two  performance  measures  relevant  to  acquisition  and  tracking  of  code 
timing  of  DS-CDMA  signals.  Since  we  are  searching  for  the  global  maximum,  we  can 
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sometimes  find  estimates  which  are  more  than  a half  chip  duration  away  from  the  true 
delay.  Such  an  estimate  is  called  an  outlier  and  is  excluded  from  the  data  before  the 
RMSE  is  calculated.  This  is  done  in  order  to  make  the  comparison  with  the  CRB 
meaningful  since  the  CRB  serves  as  a lower  bound  on  the  RMSE.  We  define  correct 
acquisition  to  be  the  event  [A  -rj  < Tc  /2 . In  order  to  measure  the  RMSE  given  correct 
acquisition,  we  define  the  RMSE  as 

RA/5£  = ^[(fi-T1)2|(ff1-T1|^rc/2)]  . (5.16) 

In  this  work,  there  is  no  tap  weight  when  a single  antenna  sensor  is  used. 

In  Figure  5.2,  the  RMSE  is  plotted  as  a function  of  the  number  of  antenna  sensors 
for  two  different  NFRs  when  M=100,  K=10  and  Eh  l /Ao=10dB  are  used.  Table  5.1 

indicates  the  relative  number  of  outlier  errors  for  the  same  situation  as  in  Figure  5.2.  It 
can  be  seen  that  the  RMSE  using  the  adaptive  antenna  array  plus  LSML  estimator  and 
MASE  algorithm  decreases  as  the  number  of  antenna  sensors  increases.  The  MASE 
algorithm  developed  to  exploit  receiver  diversity  on  flat  fading  channels  is  shown  to 
improve  the  performance  as  compared  to  that  of  LSML  timing  estimator  on  an  AWGN 
channel  where  there  exists  a DOA-dependent  phase  shift  between  two  adjacent  antenna 
elements.  In  other  words,  the  MASE  algorithm  provides  good  performance  on  a channel 
including  DOA  information  for  antenna  arrays  without  receiver  diversity  even  if  it  was 
originally  developed  for  a multiple  antenna  array  where  each  antenna  experiences 
independent  fading.  It  is  shown  that  the  pre-combining  CRB  and  post-combining  CRB 
decrease  as  the  number  of  antenna  elements  increases.  We  note  that  in  the  case  of  a single 
antenna  without  any  beamforming  process  the  CRB  derived  for  the  data  model  (3.21) 
does  not  reach  the  CRB  derived  for  the  data  model  (4.14).  In  other  words,  the  first  one  is 
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a little  larger  than  the  second  one  since  we  have  not  used  the  full  structure  of  the  data 
model  (4.14)  for  formulating  the  LSML  timing  estimator.  Another  thing  we  noticed  is 
that  the  difference  between  the  post-combining  CRB  for  OdB  NFR  and  the  pre-combining 
CRB  is  small  and  their  difference  in  the  case  of  a large  number  of  antennas  is  a little 
smaller  than  one  in  the  case  of  a small  number  of  antennas.  This  is  significant  because  it 
tells  us  that  the  loss  of  information  about  the  signal  timing  due  to  the  process  of 
beamforming  for  0 dB  NFR  is  not  much  and  there  may  also  exist  a good  timing 
estimation  method  based  on  the  beamformer  output  data. 

In  order  to  see  the  near-far  resistance  of  the  adaptive  antenna  array  plus  LSML 
timing  estimator  and  MASE  algorithm,  we  examine  the  performance  as  a function  of 
NFR  in  Figure  5.3  and  Table  5.2  when  P=3,  M=100,  K=10  and  Eb  l / Ao=10dB  are  used. 

From  Figure  5.3  which  includes  the  RMSE  for  the  case  of  a single  antenna  in  order  to 
make  a comparison,  we  can  observe  the  degradation  of  the  RMSE  for  the  high  NFR  when 
Method  1 is  used.  For  the  case  of  high  NFR,  the  use  of  an  adaptive  antenna  array 
increases  the  RMSE  of  time  delay  estimates  of  the  desired  user  and  causes  more  outlier 
errors.  Harsh  NFR  environments  could  make  the  estimated  tap  weights  poor  and  hence 
attenuate  the  desired  signal.  It  is  seen  that  the  pre-combining  CRB  is  near-far  resistant. 
Table  5.2  shows  that  when  DMI  (method  1)  is  used,  the  timing  estimation  procedure  is 
fairly  robust  to  the  near-far  problem  up  to  a NFR  of  about  12  dB.  From  this  observation, 
if  we  have  a situation  with  less  than  12  dB  NFR,  the  adaptive  antenna  array  with  3 
sensors  could  improve  the  performance  and  furthermore,  the  performance  using  Method 
1 (DMI)  is  almost  the  same  as  that  when  using  Method  2 (LCMV).  Hence  in  cases  of 
moderate  NFR,  the  additional  side  information  of  the  DOA  of  the  desired  user  does  not 
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significantly  improve  performance.  Additional  computations  of  these  methods  beyond 
that  of  the  single  antenna  sensor  based  LSML  timing  estimator  are  needed  to  calculate  the 
tap  weights  of  an  adaptive  antenna  array.  The  computations  for  estimating  the  tap  weight 
vector  of  the  adaptive  antenna  array  using  Method  1 is  mainly  due  to  the  inverse 

operation  of  a PxP  matrix,  , and  N searches  for  obtaining  the  timing  phase  of  the 
training  sequence  which  makes  the  MSE  smallest.  It  therefore  is  pointed  out  that  when 
DMI  is  used,  the  average  number  of  MATLAB  flops  for  computing  the  tap  weight  vector 
is  much  smaller  than  that  of  the  LSML  timing  estimator  if  we  use  a small  number  of 
antenna  sensors.  On  the  other  hand,  the  MASE  algorithm  requires  much  larger 
computations  than  those  of  the  LSML  timing  estimator  even  for  a small  number  of 
antenna  sensors.  For  example,  the  average  number  of  MATLAB  flops  required  by  MASE 
algorithm  with  P=2,  3 and  4 is  about  1.46,  1.92,  and  2.37  times,  respectively,  as  much  as 
that  required  by  the  LSML  timing  estimator  (P=l).  So,  the  average  number  of  MATLAB 
flops  for  calculating  the  tap  weight  vector  using  the  DMI  is  much  smaller  than  the  MASE 
requires  by  increasing  from  a single  antenna  sensor  to  multiple  antenna  sensors  and  hence 
the  computational  complexity  of  the  MASE  algorithm  is  more  expensive  than  that  of  our 
approach  using  Method  1.  It  therefore  is  shown  that  our  approach  can  improve  the 
performance  of  the  LSML  timing  estimator  with  a smaller  computational  complexity 
relative  to  that  of  the  MASE  algorithm. 

In  Figure  5.4,  the  RMSE  and  CRB  are  plotted  as  a function  of  the  number  of  users 
given  the  different  number  of  antenna  array  sensors  when  M=50,  NFR=0dB  and 
Eb, i / Aro=10dB  are  used.  Table  5.3  indicates  the  relative  number  of  outlier  errors  for  the 
same  situation  as  the  Figure  5.4.  It  is  shown  that  the  adaptive  antenna  array  offers  better 
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performance  when  a large  number  of  users  is  loaded  in  a 0 dB  NFR.  It  is  also  seen  that 
the  adaptive  antenna  array  with  3 sensors  can  increase  more  than  3 times  the  capacity 
relative  to  a single  antenna  case  with  almost  the  same  RMSE  or  CRB  performance.  It  is 
seen  that  the  pre-combining  CRB  is  almost  independent  of  the  number  of  users. 


Figure  5.2  : RMSE  using  LSML  estimator  with  random  DOAs,  M=100,  K=10  and 

£A>1/tfo=10dB 
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Table  5.1  : relative  number  of  outlier  error  with  random  DOAs,  M=100,  K=10  and 

EbxIN=\Q6B 
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Figure  5.3  : RMSE  (*)  using  LSML  estimator  with  3 antenna  sensors  (P=3),  random 
DOAs,  M=100,  K=10  and  Ebl  /Ar,=10dB 
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Figure  5.4  : RMSE  using  LSML  estimator  with  random  DO  As,  M=50,  NFR=0dB  and 
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5.1.4  Conclusions 

We  have  presented  a technique  using  an  adaptive  antenna  array  for  propagation 
delay  estimation  in  DS-CDMA  systems.  The  adaptive  antenna  array  can  suppress  the 
MAI  by  adaptively  controlling  antenna  beam  pattern  using  the  known  DOA  of  the 
desired  signal  or  the  known  training  sequence  even  if  the  desired  DS-CDMA  signal’s 
DOA  is  unknown  and  hence  the  adaptive  beamformer  output  has  less  interference.  It  is 
found  that  the  adaptive  antenna  array  helps  the  LSML  estimator  to  obtain  improved 
performance  as  compared  to  a single  antenna  based  LSML  estimator  and  the  performance 
of  the  adaptive  antenna  array  plus  LSML  timing  estimator  approaches  the  CRB  under  the 
tight  power  control  condition.  We  can  see  that  the  RMSE  and  probability  of  correct 
acquisition  of  the  adaptive  antenna  array  plus  LSML  timing  estimator  are  related  with  the 
SINR  obtained  in  chapter  3.  Considering  the  results  in  chapter  3 that  there  is  SINR 
improvement  for  0 dB  and  10  dB  NFR  whereas  20  dB  NFR  often  causes  degraded  SINR, 
the  simulation  results  in  this  Chapter  are  the  expected  ones.  That  is,  in  a less  than  about 
14  dB  NFR  situation,  the  RMSE  is  smaller  than  that  of  a single  antenna  based  LSML 
timing  estimator.  We  have  also  seen  that  the  training  sequence  based  beamforming 
(method  1)  and  DOA  based  beamforming  (method  2)  provide  almost  the  same 
improvement  under  the  tight  power  control  condition.  The  CRB  is  used  as  a performance 
bound  for  RMSE.  It  is  demonstrated  through  simulations  that  when  a large  number  of 
users  are  present  in  a perfect  power  control  condition,  our  approach  improves  the 
performance.  The  performance  of  our  approach  depends  on  an  adaptive  beamforming 
process.  The  fact  that  the  difference  between  the  post-combining  CRB  for  the  tight  power 
control  condition  and  the  pre-combining  CRB  is  small  tells  us  that  the  adaptive  antenna 
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array  followed  by  the  LSML  timing  estimator  is  a good  timing  estimation  method  under 
the  tight  power  control  condition.  It  is  shown  that  the  MASE  algorithm  gives  a little 
better  timing  estimates  at  the  cost  of  computational  complexity. 

The  existing  CDMA  standard,  IS-95  is  using  power  control  techniques  and  it  might 
be  expected  that  the  next  generation  standard  is  still  adopting  the  power  control 
techniques  due  to  their  various  benefits.  Therefore,  our  approach  may  be  applied  in 
practice. 


5.2  Timing  Estimation  for  a Flat  Fading  Channel 

In  this  section,  the  code  timing  estimation  problem  for  a single-path  ( Rk  = l)  flat- 
fading channel  is  investigated.  This  channel’s  amplitude  and  phase  are  allowed  to  vary  in 
time  even  if  the  propagation  delays  are  fixed  during  the  observation  time  interval. 

5.2.1  LSML  Timing  Estimator 

In  order  to  derive  the  LSML  timing  estimator  in  DS-CDMA  signals  over  fading 
channels,  we  assume  that  the  fading  process  is  constant  over  the  observation  time  interval 
t e [0,  MT).  Note  that  this  assumption  is  for  the  purpose  of  the  derivation  of  the  LSML 
timing  estimator. 

Then,  the  diversity  combiner  output  data  model  of  (2.58)  can  be  written  as 

y(m)  = Dj  Zj  (m)  + e(m) 

=iuA,(Ti)Zi(ra)+e(m)6Cw  , (5-17) 

where  Dj  =xuA1(r1)  with  jcu  = V^fwfQu. 
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Now,  we  can  derive  the  LSML  timing  estimator  of  (5.1)  from  the  data  model  of 
(5.17)  where  the  contribution  of  all  undesired  signals  is  lumped  in  with  that  of  the 
thermal  noise. 

5.2.2  MASE  Algorithm 

In  order  to  develop  the  MASE  algorithm  [26],  we  need  the  same  assumption  as  that 
for  deriving  the  LSML  timing  estimator  for  fading  channels.  That  is,  the  fading  process  is 
assumed  to  be  constant  over  the  observation  time  interval  te  [0,  MT)  even  if  the 
amplitude  and  phase  of  the  received  signal  is  time-varying.  Remember  that  this 
assumption  is  for  the  purpose  of  developing  the  MASE  algorithm. 

Note  that  the  MASE  algorithm  is  based  on  the  NPxl  array  output  vector  as  in 
(4.1 1)  from  P antenna  sensors  whereas  the  LSML  timing  estimator  employing  a diversity 
combiner  is  based  on  the  diversity  combiner  output  vector  of  (5.17).  The  MASE 
algorithm  is  given  by  (5.9). 

5.2.3  Numerical  Results 

The  performance  of  the  estimator  in  each  simulation  was  obtained  from  1000 
Monte-Carlo  runs.  The  fading  process  is  generated  by  passing  white  gaussian  noise 
through  a third  order  low  pass  filter.  The  data  rate  is  assumed  to  be  9600  bps.  The 
performance  of  a diversity  combiner  based  LSML  timing  estimator  on  a time- varying 
fading  channel  is  evaluated  by  assuming  single  path  (Rk  = l)  flat  fading.  We  consider 

only  Method  1,  based  on  training  sequences  to  calculate  the  tap  weights  of  the  adaptive 
antenna  array  on  the  time- varying  fading  channel  because  the  DOA  based  beamforming 
may  not  be  effective  due  to  the  time- varying  phase  of  the  fading  process.  This  adaptive 


antenna  array  is  referred  to  as  a diversity  combiner.  We  also  compare  the  diversity 
combiner  based  LSML  timing  estimator  with  the  MASE  algorithm  in  [26]. 
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In  Figure  5.5,  the  RMSE  of  the  diversity  combiner  based  LSML  timing  estimator 
and  the  MASE  algorithm  is  plotted  as  a function  of  the  number  of  antennas  when  M=50, 
K=20,  NFR=0  dB,  Eb  l / Ao=10dB  and  fd  (fading  rate)  = 70  Hz  are  used.  Table  5.4 

indicates  the  relative  number  of  outlier  errors  for  the  same  simulation  condition  as  Figure 
5.5.  It  is  shown  that  the  RMSE  and  the  relative  number  of  outlier  errors  of  both 
approaches  decrease  as  the  number  of  antenna  sensors  increases.  It  is  seen  that  the  MASE 
algorithm  gives  a little  better  timing  estimates  at  the  cost  of  a computational  complexity. 
The  use  of  multiple  antenna  sensors  increases  the  probability  of  correct  acquisition  for 
the  LSML  timing  estimator. 

In  order  to  examine  the  sensitivity  of  the  diversity  combiner  plus  LSML  timing 
estimator  and  MASE  algorithm  to  the  fading  rate,  we  plot  the  performance  as  a function 
of  the  fading  rate  in  Figure  5.6  and  Table  5.5  when  M=50,  P=4,  K=20,  NFR=0  dB  and 
£*1  /jVo=10dB  are  used.  Even  if  the  performance  of  RMSE  and  probability  of  correct 

acquisition  for  both  approaches  is  degraded  as  the  fading  rate  increases,  the  performance 
of  correct  acquisition  is  still  good  for  the  given  scenario. 

5.2.4  Conclusions 

The  comparison  of  the  performance  of  the  diversity  combiner  plus  LSML  timing 
estimator  with  that  of  the  MASE  algorithm  has  been  studied  on  a single-path  flat-fading 
channel.  The  diversity  combiner  can  suppress  the  interfering  signals  by  exploiting  the 
spatial  diversity  of  the  received  signals.  Thus,  the  diversity  combiner  based  LSML  timing 
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estimator  can  improve  the  performance  of  LSML  timing  estimator.  It  is  shown  that  the 
MASE  algorithm  gives  a little  better  timing  estimates  at  the  cost  of  computational 
complexity. 


Figure  5.5  : RMSE  on  a time-varying  flat  fading  channel  with  M=50,  K=20,  NFR=0dB, 

Eb  i / Ao=10dB  and  fd=70Hz 


P=1 

P=2 

P=3 

P-4 

P=6 

P=10 

Diversity  combiner  based 
LSML 

0.066 

0.011 

0.001 

0 

0 

0 

MASE 

0.066 

0.003 

0 

0 

0 

0 

Table  5.4  : relative  number  of  outlier  error  on  a time- varying  flat  fading  channel  with 
M=50,  K=20,  NFR=0dB,  Ebl  / Ao=10dB  and  fd=70Hz 
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Figure  5.6  : RMSE  on  a time- varying  flat  fading  channel  with  M=50,  K=20, 

Eb  i / Aro=10dB  and  NFR=0dB 


Fading  rate  (Hz) 

10 

40 

70 

100 

130 

160 

190 

P=1  ( LSML  ) 

0.010 

0.027 

0.066 

0.138 

0.212 

0.310 

0.411 

P=4 

Diversity 
combiner  based 
LSML 

0 

0 

0 

0.001 

0.001 

0.010 

0.031 

MASE 

0 

0 

0 

0 

0 

0 

0 

Table  V : relative  number  of  outlier  error  on  a time- varying  flat  fading  channel  with 
M=50,  K=20,  Eb  ] / A^o=10dB  and  NFR=0dB 
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5.3  Timing  Estimation  for  a Frequency  Selective  Fading  Channel 

In  the  previous  section  we  have  seen  that  both  the  MASE  algorithm  and  the 
diversity  combiner  plus  the  LSML  timing  estimator  work  well  for  a single  path  time- 
varying  flat  fading  channel  under  a perfect  power  control  condition.  In  this  section,  we 
propose  an  efficient  algorithm  for  an  approximate  maximum  likelihood  approach  of 
jointly  estimating  the  code  timings  of  a desired  user  for  DS-CDMA  systems  that  consist 
of  an  arbitrary  antenna  array  at  the  receiver  and  works  in  the  presence  of  frequency 
selective  fading  and  the  near-far  problem.  The  algorithm  is  referred  to  as  the  extended 
MASE  (EMASE).  It  is  shown  that  the  EMASE  provides  good  performance  in  frequency 
selective  fading  channels  and  in  the  presence  of  the  near-far  problem.  The  acquisition 
performance  of  the  MUSIC  timing  estimator  with  a single  antenna  is  not  good  for  small 
observation  windows  [47]  and  also  low  signal  to  noise  ratios  (SNR)  [22].  A possible 
reason  for  the  poor  acquisition  performance  for  small  observation  windows  is  that  the 
path,  whose  timing  is  to  be  estimated,  might  be  in  a deep  fade  during  the  observation 
interval  which  results  in  a very  noisy  estimate.  In  a stationary  channel,  much  larger 
observation  intervals  could  be  used  to  improve  the  estimates  at  the  cost  of  additional 
complexity.  When  the  channel  changes  dramatically,  the  interval  can  not  be  increased 
arbitrarily  because  the  system  may  need  to  adapt  to  time-varying  fading  channels  more 
quickly.  Thus,  by  exploiting  spatial  diversity  via  multiple  antennas  its  acquisition 
performance  even  for  small  observation  windows  or  low  SNR  could  be  improved  as 
compared  to  that  of  the  MUSIC  with  a single  antenna.  So,  the  MUSIC  timing  estimator  in 
[47]  is  extended  to  the  case  of  using  a multiple  antenna  array  for  estimating  multipath 
timings  in  frequency  selective  fading  channels.  Our  approach  is  different  from  one  used 
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in  [25]  in  that  we  use  the  spatial  diversity  to  allow  for  rapid  acquisition  while  the  authors 
of  [25]  use  the  spatial  diversity  to  achieve  higher  capacity.  The  multiple  antennas  based 
MUSIC  timing  estimator  is  referred  to  as  MAMUSIC.  It  is  shown  that  MAMUSIC 
performs  significantly  better  than  the  single  antenna  based  MUSIC  timing  estimator, 
especially  when  the  number  of  observation  bits  is  small  or  SNR  is  low.  The  performance 
of  the  MAMUSIC  estimator  is  compared  to  that  of  the  EMASE  via  simulation  results. 
These  extended  timing  estimation  algorithms  for  frequency  selective  fading  channels  are 
presented  in  the  following  subsection. 

5.3.1  Timing  Estimators 

We  develop  two  multipath  timing  estimation  methods  using  spatial  diversity  via 
multiple  antennas.  The  first  one  is  the  EMASE  algorithm  which  has  been  modified  for 
frequency  selective  fading  channels  and  the  second  one  is  the  MAMUSIC  timing 
estimator  which  now  exploits  spatial  diversity. 

5.3. 1.1  EMASE  algorithm 

We  assume  that  the  fading  process  is  constant  over  the  observation  time  interval 
t e [0,  MT ) . Note  that  this  assumption  is  for  the  purpose  of  the  derivation  of  the  timing 
estimator.  Hence,  the  derived  EMASE  estimator  is  an  approximate  ML  estimator  for  a 
time-invariant  system,  however,  simulation  experiments  are  conducted  for  time-varying 
fading  scenarios  to  examine  sensitivity  to  the  fading  rate.  The  EMASE  algorithm  is 
derived  by  using  a known  training  signal  and  modeling  the  multiple  access  interference 
(MAI)  and  the  additive  noise  as  an  unknown  colored  Gaussian  random  process.  Thus,  the 
received  vector,  r (p\m)  of  (2.39)  can  be  expressed  as 
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where 


r(p)  (m)  = (d^Zj  (m)  + e(p)  ( m))T , 


®iP)  _ Ai(riti,Ti2»,"»T1B  )b|^ 


(5.18) 


(5.19) 


and  the  vector  e(p)(m)  is  assumed  to  be  a complex  Gaussian  random  vector  with  zero- 


mean  and  arbitrary  covariance  matrix  Q with 


= Q , 


(5.20) 


By  defining  r(m)  as  in  (4. 1 1),  the  joint  probability  distribution  of  the  received 
vectors  is  then  given  by 


/[{?(«)}",!  | D^.q] 

( PM 


= * ,PU  exP^-XZ|' r(p)T(m)-D[p)z !(m)l  Q T r^^-D^z ,(/«)] 

71  | Q|  [ p=\  m=l L J L J 


(5.21) 


We  can  then  derive  a joint  timing  estimation  procedure  in  a manner  similar  to  what  was 

done  for  the  MASE  [26].  We  form  the  unstructured  estimates 

d''>=r“  R;'=[dj'>  d<'>] 

(5.22) 

/>=i 

(5.23) 

where 

R*,  =-^Zzi  ^yp)\m) 

M m=l 

(5.24) 

Rr,r,  =~-^r(p)T  (m)r(p)\m) . 

M m= 1 

(5.25) 

Then,  we  seek  to  minimize  the  cost  function 
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C,  = 


-A;„A-r  Q_i(dfp) 


(5.26) 


The  value  of  which  minimizes  this  cost  function  is 


A"’  =(a;J'q-,a|j 

Then,  the  estimates  of  r = { r, , , 
cost  function 


+ a;"q-1a;i2 

ri,2 » " ■ » ri,/?,  } 


Y (x'“Q-ld\p)  + A'^Q-’d^ ).  (5.27) 

will  be  those  values  which  minimize  the 


C2  - Cj 


(5.28) 


The  equation  in  (5.28)  can  be  shown  to  be  equivalent  to  maximizing  the  cost  function 


C3 


+ A'»Q-'d</> 


)1 


p= 1 


(5.29) 


This  joint  estimate  requires  an  -dimensional  optimization  over  t and,  hence, 
expensive  computations.  Therefore,  we  present  a simple  algorithm  for  trying  to  maximize 
the  cost  function.  The  approach  is  to  first  optimize  the  one-dimensional  cost  function  to 
obtain  one  path’s  time  delay  (perhaps,  the  strongest  one).  Then,  the  maximization  of  two- 
dimensional  cost  function  is  performed  to  get  the  second  path’s  timing  with  the  first 
dimension  fixed.  Then,  we  continue  in  this  manner  for  as  many  paths  as  desired.  Let 
C3 A (r, , , r,  2 , • • • , r,  Rt ) be  the  cost  function  defined  previously  based  on  an  -path 
channel  model.  We  can  obtain  f,  ,,  f12, •••,  f1Si  by  solving  the  following  R{  one- 
dimensional maximization  problems.  That  is, 


f,,  = arg  max  C3l  (r) 


(5.30) 
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and  then 


fj  2 = arg  max  C32  (f  u ,t)  (5.3 1) 

T 

and  then,  finally 

rlA  = arg  max  C3^  *”^"1.^—1  »t).  (5.32) 

Now,  given  the  first  one-dimensional  search  of  C3l  (r)  we  develop  two  different 
procedures  to  maximize  the  next  cost  function  C3  2 (fu  ,t)  and  then 
Q,r  (fu , f i,2 . • • • , Ti,r-i , f ) where  r = 3, 4,  --•,/?! . 

5.3.1. 1. 1 Method  1.  Given  the  first  one-dimensional  search  of  C3 4 (t)  which  can  be 

accomplished  by  rooting  a second-order  polynomial  we  could  consider  a procedure  to 
maximize  the  next  cost  function  C3  2(tu,t).  Let  ru  = puTc  + pnTc  and  r = pTc  + pTc . 

We  can  get  the  C32( ru,r)  after  an  analytical  matrix  calculation,  as  described  in 
Appendix  B. 


C^X)~~5^) 

+«i('Vw  +niP)p+n(3p)pilp+n(/)pl1  + n\p)p2  +n(6p)pl1p+n(/)pup2  +nlP)p2lp2) 


(5.33) 


^0+^1  A*u  +d2P+d3  pup  + d4pl1+dip2+  d6  p^p + dn  pup 2 +dtpl,p2 


where  the  coefficients  n\p)  and  dt  for  i = 0,l,---,8  are  given  in  Appendix  B. 

The  C3  2(tu,t)  is  in  the  form  of  a ratio  of  second-order  polynomials  in  /i, , and  p where 
Mu  =SU /Tc  e [0,1 ) and  p = S/Tc  e [0,1 ).  Since  ru  can  be  replaced  by  fu , the 
C3  2(tii,t)  can  be  written  as 


m0  +mlp  + m2p 2 
32  U’  oifiu’L1)  b0+bxp+b2p 2 


(5.34) 
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where 


=S(/Io/,)  +n[P)Pu  +«4P>Ai)» 

p= 1 

bo  =d0  + d{  fin  + d+fin 

(5.35) 

m,=X(<>+<>A4+»rA;1). 

P= 1 

b\  — d2  + d3  filtl  + d6  fin 

(5.36) 

m2  = X + niP)P  1.1  + ^ fit  )’ 

P= 1 

bi  -d5  + dj  fin  + ds  fiu  • 

(5.37) 

Thus,  the  time  delay  of  the  second  path  can  be  obtained  by  also  rooting  a second-order 
polynomial  in  /r  of  (5.34).  Then,  we  can  continue  in  this  way  for  the  next  path  by 
obtaining  the  cost  function  C3  r(Tu,T12,"-,Ti>r_1,T)  with  a ratio  of  polynomials  in  jiu, 
/r12,-”  and  /rlr  where  fil  r =Sl  r/Tc  e [ 0, 1 ) . However,  in  general  it  seems  too  complex  to 
obtain  the  cost  function  C3  i.(t11,t12,-*-,t1i,._1,t)  with  a ratio  of  polynomials  in  /ru, 

/r12,---  and  nu  because  it  involves  the  analytical  matrix  inverse  computation  of  the  rxr 
matrix  A^Q_1A^  + A|"Q_1Aj  2 . That’s  why  we  consider  only  two-path  timing 

estimation  for  this  method  in  the  simulations  even  if  the  system  model  is  based  on  a three 
path  fading  channel.  Therefore,  in  the  next  subsection  we  present  another  method  to 
easily  extend  to  the  case  of  more  paths. 

5.3. 1. 1 .2  Method  2.  We  present  the  simpler  procedure  to  optimize  the  cost  functions 

Based  on  the  r-path  channel  model,  we  follow  a similar 
procedure  as  the  one  taken  from  equation  (5.26)  to  equation  (5.29).  Given  the 
Wu 311(1  PiP) > $$ .....  A(,r-i > the  Q in  (5.26)  can  be  written  as 
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p- 1 


i=l 


( r— 1 


i=i 


»-i 


d,(P>  - XaU  feu  JAP  ~ 8l,  fel, 

1=1 


where 


a(p)  ^u(fij)^p)+aiH2(rij)d(2p) 

au  feu  K fe  u )+  ai  2 feu  )ali2  feu ) 


■u  feu  )=  Q"X Iau  (PuTc ) ai,i  (tPu  + 1FC )] 


1a 

-Pu 

A; 


, /=l,2,-,r-l 


d,^  =Q"^d|'\  t = l»  2 . 
Then,  the  cost  function  can  be  written  as 


2 

r , / r-l  > 

1 

2 

I 

ssfcj  a,<«  -£*;;> 

i=l 

L l '=i  J 

J 

'=1  ||al,l  fel.r  ]j|  +|aufeul' 


where 


V?  ~ au  feu  W ■ 


We  can  then  rewrite  C3_,  in  (5.43)  as 


P H/x 


where 


(5.38) 

(5.39) 

(5.40) 

(5.41) 

(5.42) 

(5.43) 

(5.44) 


P = [ 1-P  p]r 


(5.45) 
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( p 

f 

f 

r-\  > 

( 

r— 1 V 

I 

ah 

Al,l 

d{rt 

-1^ 

+ A" 

d(2p) 

l'-' 

/-I  ) 

/=i  7. 

r 

( r- 1 > 

( r-1  V 

a£ 

, 7H 
+ Ai,2 

a«'>-£v« 

- 

k /= i 7 

k /=i  7. 

7 

(5.46) 


H = Re{A"Au+A"A1>2}  (5.47) 

Aw=Q'X[a1,(prc)  au([p  + l]r.)],  / = 1,2.  (5.48) 

So,  the  problem  becomes  rooting  a second-order  polynomial. 

5.3. 1.2  MAMUSIC  timing  estimator 

The  EMASE  algorithm  requires  knowledge  of  a training  sequence.  In  this 
subsection,  the  MUSIC  estimator,  [47]  which  does  not  need  training,  is  extended  to  the 
case  of  using  a multiple  antenna  array  for  estimating  multipath  timings  in  frequency 

selective  fading  channels.  The  correlation  matrix,  R,  for  r(p\m)  is 


R = E 


= ASA"  +ct2In, 


(5.49) 


where  S = £ s(/’)(7n)s(p)//(w) 


. We  assume  that  A and  S have  full  rank  and  the  signal 


subspace  and  noise  subspace  are  defined  as  the  column  space  of  A,  i.e.,  range( A),  and  the 
orthogonal  complement  to  the  signal  space,  respectively.  Furthermore,  since  A has  rank 
2 R,  the  signal  subspace  will  have  dimensionality  2 R.  We  can  form  an  eigen- 
decomposition  of  R as 


R = E!A1Ef  +cr2E„Ef , (5.50) 

where  As  is  a diagonal  matrix  of  2 R largest  eigenvalues  of  R and  Es  is  a matrix  of  the 
eigenvectors  of  R corresponding  to  the  2 R largest  eigenvalues  of  R and  E„  is  a matrix  of 
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the  eigenvectors  of  R corresponding  to  the  N-2R  smallest  eigenvalues.  The  columns  of 
Ej  span  the  signal  subspace  while  the  columns  of  E„  span  the  noise  subspace. 

Since  R is  unknown  to  the  receiver,  it  is  estimated  by  the  sample  correlation  matrix 


=-^Vy’rw’'(m)rl'’>'(m)  = E,A,Ef  +<t2E.E,, 


PM 

p= l m=i 


(5.51) 


where  the  columns  of  E„  are  the  eigenvectors  corresponding  to  the  N-2R  smallest 
eigenvalues  of  RPM  . Note  that  when  P=l,  MAMUSIC  is  reduced  to  MUSIC.  In  an 
alternative  approach  used  in  [25]  the  correlation  matrix  is  defined  as  is[r(/w)rH(ffj)] 

1 M 

where  the  sample  autocorrelation  matrix  is  computed  by  — V r(m)r//  ( m ) . This  allows 

M 1 

for  synchronizing  to  more  users  than  a single  antenna  based  MUSIC  estimator  by 
increasing  the  dimensionality  of  the  received  vector  due  to  the  presence  of  multiple 
antennas. 

Since  the  known  desired  signal  is  contained  in  the  signal  subspace,  it  must  be 
orthogonal  to  the  noise  subspace  and  hence  the  timing  estimator  is  to  be  the  value  of  the 
timing  for  which  the  known  code  sequence  is  nearest  to  being  orthogonal  to  an  estimate 
of  the  noise  subspace.  Let  the  MAMUSIC  cost  function  for  the  first  user  be  defined  as 


J,M= 


■SM*. 

2 

+ 

KWE.f 

II  H ( ^ 

I2 

H ( "Jl2 

|*u(u| 

+ 

|aU  (U|| 

(5.52) 


The  MUSIC  cost  function  in  (5.52)  is  slightly  different  from  the  one  used  in  the  original 
work  [20].  It  results  in  a computationally  simpler  algorithm  even  if  there  may  be  some 
degradation  in  using  the  modified  cost  function.  It  is  easy  to  show  that  J,  (r)  is  a rational 
function  of  two  second-order  polynomials  in  fi  for  each  chip  interval.  Let  3 denote  the 
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set  of  costs  corresponding  to  the  candidate  estimate  of  t . We  now  describe  the  procedure 
for  finding  multiple  timing  estimates  from  the  MUSIC  algorithm  as  the  following  steps: 
Stepl:  Let  3 = {j1(0)> -1FC)} 

Step  2:  For  p = 0,l,---,N 

(a)  Find  the  roots  fLx  and  p2  of  a second-order  polynomial. 

(b)  If  a root  /2,  e [0,  1 ),  then  add  Ji  {[p  + ft,  Jr,. ) to  3 . 

Step  3 : The  timing  estimates,  ^ } , are  found  as  the  members  of  the 

set  of  candidate  timing  estimates  corresponding  to  the  R{  smallest  values  of 
the  cost  function. 

5.3.2  Numerical  Results 

The  simulated  system  uses  Gold  sequences  with  N=31  chips  per  bit  operating  over 
a fading  channel.  The  performance  of  the  estimator  in  each  simulation  was  obtained  from 
1000  Monte-Carlo  runs.  The  paths  from  each  user  are  assumed  to  have  the  same  mean 
power.  The  timing  offset  zk  r and  data  bits  of  all  users  are  independent  of  each  other.  The 

fading  process  is  generated  by  passing  white  Gaussian  noise  through  a third  order  low 
pass  filter.  In  these  simulations,  the  data  rate  is  assumed  to  be  9600  bps.  By  assuming  a 
three-path  frequency  selective  fading  model,  the  performance  of  the  EMASE  and 
MAMUSIC  on  a time- varying  fading  channel  is  evaluated  in  two  cases.  The  first  one  is 
that  if  at  least  one  of  fu , f1>2  and  f13  produced  by  the  timing  estimator  is  greater  than  a 

half  chip  duration  away  from  the  true  delay  of  the  corresponding  path,  there  is  an 
incorrect  acquisition.  This  case  is  referred  to  as  3 -path  in  the  following  plots.  The  second 
one  is  that  we  obtain  the  timing  estimates  of  two  paths  (perhaps,  two  strongest  paths) 
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among  three  paths  of  the  desired  user’s  signal  even  if  the  system  model  is  based  on  three 
paths  fading  channel  This  case  is  referred  to  as  2-path  in  the  following  plots  where  there 
is  a correct  acquisition  if  all  the  timing  estimates  of  two  paths  are  less  than  a half  chip 
duration  away  from  the  true  delay  of  the  corresponding  path.  The  propagation  delays  of 
each  user’s  signal  are  randomly  chosen  and  then  fixed  to  be 
fru  t12  /Tc  =[8.3417  11.7818  16.51 19f . 

Figures  5.7  and  5.9  show  the  probability  of  correct  acquisition  of  EMASE  and 
MAMUSIC  as  a function  of  the  NFR  with  M=50,  K=4,  Ebx  /No=10dB  (each  path  of  the 
desired  user’s  signal)  and  fd=100Hz  (fading  rate)  for  three  different  numbers  of  antennas. 
Figure  5.8  shows  the  root  mean  square  error  (RMSE)  of  the  first  path  given  a correct 
acquisition  in  the  same  situation  as  in  Figure  5.7.  The  performance  is  shown  to  be 
significantly  improved  as  the  number  of  antennas  increases.  It  is  seen  that  the  EMASE 
provides  better  acquisition  performance  as  compared  to  that  of  the  MAMUSIC.  Note  that 
the  difference  of  an  acquisition  performance  between  the  EMASE  and  MAMUSIC  for  1 
and  2 antenna  cases  is  significant  whereas  the  acquisition  probability  of  the  MAMUSIC 
for  4 antennas  is  greater  than  90  %.  The  EMASE  and  MAMUSIC  are  shown  to  be  fairly 
robust  to  the  near-far  problem.  Note  that  the  CRB  is  near-far  resistant. 

In  Figures  5.10,  5.11,  and  5.12,  the  performance  is  plotted  as  a function  of  the 
number  of  observation  bits  when  Eb  l / Ao=10dB  (each  path),  K=4,  NFR=10dB  and 

fd=100Hz  are  used.  The  use  of  larger  observation  bits  does  not  improve  the  acquisition 
performance  of  EMASE.  The  assumption  that  the  fading  process  is  constant  over  the 
observation  time  interval  was  used  for  the  derivation  of  the  EMASE  algorithm  which  is 
based  on  the  estimates  of  D,(rt  from  (5.22)  and  of  (5.39)  which  are  a function  of  the 
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fading  process.  Because  in  fact  the  fading  process  is  time  varying,  these  fixed  estimates 
can  not  accurately  represent  the  time- varying  quantities  as  the  length  of  the  observation 
window  increases.  That’s  why  the  performance  of  the  EMASE  does  not  improve  as  the 
number  of  observation  bits  increases.  On  the  other  hand,  the  MAI  and  noise  covariance 
matrix  is  estimated  as  the  average  of  those  estimated  from  each  antenna  output,  as  shown 
in  (5.23).  So,  the  performance  of  the  EMASE  is  improved  by  using  the  multiple  antennas. 
It  is  shown  that  the  performance  of  the  MAMUSIC  estimator  improves  significantly  as 
the  number  of  observation  bits  increases.  Note  that  the  performance  improvement  is 
significant  by  using  multiple  antennas  even  for  a small  number  of  observation  bits.  We 
need  about  60  observation  bits  for  4 antennas  in  order  to  get  about  95  % correct 
acquisition  probability  for  the  MAMUSIC  estimator  for  3 -path  which  requires  on  the 
order  of  PMN2  (=  240 N2 ) multiplications  to  compute  a sample  correlation  matrix.  This 
same  acquisition  performance  needs  about  280  observation  bits  for  the  MUSIC  estimator 
with  a single  antenna,  which  requires  on  the  order  of  MN2  (=  280N2 ) multiplications. 
Therefore,  in  order  to  obtain  almost  the  same  performance,  the  computation  of  the 
MUSIC  estimator  can  be  slightly  reduced  by  using  MAMUSIC.  Note  that  the 
MAMUSIC  estimator  provides  better  performance  than  EMASE  for  3 -path  when  more 
than  120  observation  bits  are  used. 

Figures  5.13  and  5.14  show  the  sensitivity  of  the  acquisition  performance  of 
EMASE  and  MAMUSIC  to  the  fading  rate  with  M=50,  K=4,  Ebx  /No=10dB  (each  path) 

and  NFR=6dB.  The  EMASE  with  4 antenna  sensors  is  shown  to  be  relatively  insensitive 
to  the  fading  rates  used  in  this  simulation.  Note  that  the  performance  of  the  MAMUSIC 
estimator  slowly  improves  as  the  fading  rate  increases.  In  a slowly  fading  environment. 
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the  path  whose  delay  is  to  be  estimated,  might  be  in  a deep  fade  during  the  50  observation 
bits  which  could  make  the  estimation  difficult.  In  the  fast  time-varying  case,  the  entire 
observation  window  will  most  likely  not  be  in  a deep  fade,  hence  a better  estimates  of  the 
noise  subspace  can  be  obtained.  Therefore,  the  MAMUSIC  estimator  is  more  appropriate 
for  rapidly  time- varying  channels. 

The  performance  of  EMASE  and  the  MAMUSIC  estimator  is  plotted  as  a function 
of  Eb  l/  Na  (each  path)  in  Figure  5.15,  5.16  and  5.17  when  M=200,  K=4,  NFR=10dB  and 
fd=100Hz  are  used.  It  is  seen  that  the  use  of  multiple  antennas  provides  much  better 
performance  than  a single  antenna  case.  The  EMASE  (P=4)  works  even  for  low  Ebl  / N„ 

and  gives  better  performance  than  the  MAMUSIC  (P=4).  When  Eb  l/N0  is  above  8 dB, 
the  performance  of  the  MAMUSIC  estimator  is  better  than  that  of  EMASE. 

5.3.3  Conclusions 

Two  different  timing  estimators  of  DS-CDMA  signals  exploiting  spatial  diversity 
are  proposed  for  estimation  of  multipath  timings  in  a frequency  selective  fading  channel 
and  near-far  environment.  The  proposed  EMASE  estimator  with  4 antenna  sensors  is 
shown  to  be  fairly  robust  to  the  near-far  problem  and  relatively  insensitive  to  the  fading 
rate.  The  EMASE  also  works  well  for  the  small  observation  windows  and  low  Ebl/  N„. 

Furthermore,  the  proposed  algorithm  for  an  approximate  ML  estimator  does  not  require  a 
multi-dimensional  optimization  and  the  code  timings  of  all  paths  of  a desired  user  are 
obtained  by  rooting  a second-order  polynomial,  which  is  computationally  efficient.  The 
MAMUSIC  has  also  near-far  resistance.  The  MAMUSIC  estimator  significantly 
improves  the  performance  relative  to  the  single  antenna  MUSIC  for  small  observation 
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times  and  low  Eb  l / Na . The  MAMUSIC  estimator  can  obtain  almost  the  same 

performance  as  that  of  the  single  antenna  MUSIC  estimator  without  additional 
computation  complexity.  Both  EMASE  and  MAMUSIC  provide  much  better  acquisition 
performance  for  2-path  estimation  than  for  3 -path  over  three  paths  time- varying  fading 
channels. 
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Figure  5.7  : Probability  of  correct  acquisition  for  3-path  using  Method  2 of  EMASE  and 
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Figure  5.8  : RMSE  of  the  time  delay  of  the  first  path  for  3-path  using  Method  2 of 

EMASE  and  MAMUSIC 
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Figure  5.9  : Probability  of  correct  acquisition  for  2-path  using  EMASE  and  MAMUSIC 
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Figure  5.10  : Probability  of  correct  acquisition  for  3-path  using  Method  2 of  EMASE  and 

MAMUSIC 


Figure  5.11:  RMSE  of  the  time  delay  of  the  first  path  for  3 -path  using  Method  2 of 

EMASE  and  MAMUSIC 
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Figure  5.12  : Probability  of  correct  acquisition  for  2-path  using  EMASE  and  MAMUSIC 
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5.13  : Probability  of  correct  acquisition  for  3 -path  using  Method  2 of  EMASE  and 

MAMUSIC 


Figure  5.14  : Probability  of  correct  acquisition  for  2-path  using  EMASE  and  MAMUSIC 
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Figure  5.15  : Probability  of  correct  acquisition  for  3-path  using  Method  2 of  EMASE  and 

MAMUSIC 


Figure  5.16  : RMSE  of  the  time  delay  of  the  first  path  for  3-path  using  Method  2 of 

EMASE  and  MAMUSIC 
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Figure  5.17  : Probability  of  correct  acquisition  for  2-path  using  EMASE  and  MAMUSIC 


CHAPTER  6 

MMSE  DETECTION  WITH  MULTIPLE  ANTENNA  SENSORS 


This  chapter  investigates  MMSE  detection  with  multiple  antenna  array  sensors.  In 
chapter  3,  it  has  been  shown  that  the  adaptive  antenna  array  or  diversity  combiner  can 
reduce  the  MAI  and  hence  improve  the  SINR  performance.  The  MMSE  detector  using 
spatial  diversity  via  multiple  antenna  sensors  may  reduce  the  effect  of  channel  fading  and 
suppress  the  MAI  and  hence  have  more  improved  BER  and  capacity  than  that  of  a 
MMSE  detector  without  a multiple  antenna  array.  In  this  chapter,  we  examine  how 
spatial  diversity  via  a multiple  antenna  array  is  used  for  MMSE  detection  and  how  much 
the  multiple  antenna  array  could  improve  the  BER  and  capacity  performance  of  the 
MMSE  detector  on  time- varying  fading  channels. 


Assuming  the  receiver  has  the  knowledge  of  the  propagation  delay  of  all  paths  of 
the  desired  signal,  rl 4 can  be  set  to  zero  without  loss  of  generality.  The  conventional 

receiver  for  DS-CDMA  would  correlate  the  received  vector  for  the  mth  data  bit  with  a 
local  replica  of  the  desired  user’s  spreading  code.  That  is,  a decision  on  the  wth  bit  in  a 
flat  fading  channel  would  be  of  the  form 


6.1  Minimum  Mean-Squared  Error  Receiver 


(6.1) 


98 


99 


Note  that  the  detector  must  track  and  account  for  the  phase  on  the  desired  user’s  signal. 
This  decision  is  optimum  for  the  single-user  case  but  can  lead  to  substantially  degraded 
performance  in  the  presence  of  multiple-access  interference. 

Figure  6. 1 depicts  the  conventional  MMSE  receiver.  The  adaptive  algorithm  adjusts 


The  error  signal,  eft  ( m ) , is  formed  as  the  difference  between  the  filter  output  and  the 

desired  user’s  bit.  The  receiver  will  therefore  require  that  a known  data  sequence,  or 
training  sequence,  is  transmitted  before  any  actual  data  can  be  transmitted.  When  the 
MSE  has  lowered  to  an  acceptable  level,  data  transmission  begins  and  the  error  signal  is 
formed  as  the  difference  between  the  filter  output  and  the  hard  decision. 


These  modes  are  referred  to  as  training  and  decision  directed  modes,  respectively. 


the  receiver  vector  which  minimizes  the  mean-squared  error  (MSE),  J = 


4h(>>|2} 


(m)ra)T  (m) 


(6.2) 


r(1)(m)  Adaptive 

^ Filter 

w 2,1  (m) 


> Re{  } ► Sgn  () 


dx(m) 


Adaptive 

algorithm 


dx  (m) 

Figure  6.1:  Block  diagram  of  a minimum  mean  square  error  receiver 
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It  is  well-known  that  the  MSE  can  be  written  as 


= ^ |gu  (m>|2  j = ^ 


di  ( m ) - w 2jH  (m)r(1)r  (m) 


(6.3) 


The  receiver  vector  which  minimizes  the  MSE  is  well  known  to  be 


w 24  (w)  = R (1)_1  (m)p  a)  (m) , 


(6.4) 


where  R(1)(/n)  = £^rfl)r(m)r(1)  (m)J  and  pa)(m)  = £^d1(m)ra)r(m)J  are  the  correlation 


matrix  and  steering  vector,  respectively.  Typically,  the  receiver  vector  can  be  updated  by 
the  least  mean  squared  (LMS)  or  the  recursive  least  squared  (RLS)  algorithm. 

The  conventional  MMSE  detectors  may  experience  the  following  problems  on 
Rayleigh  fading  channels  [53] : (i)  locking  180  degrees  out  of  phase  to  the  correct  user’s 
signal ; and  (ii)  false  locking  to  any  of  the  interferer’s  signals,  which  can  occur  during  a 
fade.  The  first  problem  can  be  solved  by  using  differential  encoding  and  detection.  The 
second  problem  makes  the  adaptive  MMSE  receiver  useless  on  a Rayleigh  fading 
channel.  In  [54],  two  different  adaptive  algorithms  for  the  MMSE  detector  in  a single- 
path flat  fading  channel  are  considered.  Another  alternative  is  the  modified  MMSE 
receiver  [53]  in  Figure  6.2  where  phase  variations  of  the  desired  user’s  signal  are 
removed  from  the  received  signal  prior  to  entering  the  adaptive  filter.  A separate  Lth 
order  linear-predictor  is  used  to  track  the  channel  (phase  and  amplitude)  prior  to  the 
adaptive  filter. 

Since  the  conventional  MMSE  receiver  loses  phase  lock  during  deep  fades,  the 
phase  precompensation  makes  the  MMSE  receiver  perform  properly  even  in  a Rayleigh- 
fading  channel.  The  modified  MMSE  detector  in  a flat  fading  channel  forms  its  data 


decisions  as 
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di(m)  = sgn 


= sgn 


^w«>r(ra)Re[e-ii”<-,r<1>!' 


where  r $)(m)  = e 


-*?»r<nr 


(m). 


(6.5) 


Figure  6.2  : Block  diagram  of  the  modified  MMSE  receiver 

Since  the  adaptive  filter  may  lock  180°  out  of  phase  to  the  correct  user’s  signal, 
differential  encoding  and  decoding  are  used  for  the  rest  of  this  Chapter,  as  shown  in 
Figure  6.3.  The  mth  differentially  encoded  data  bit  for  the  kth  user  can  be  written  as 

d'k(m)  = dk(m)d'l(m-l),  (6.6) 

where  dk  (m)  is  the  actual  data  bit.  For  BPSK  signal,  the  decision  on  the  mth  data  bit  is 
made  according  to 

rf,  (m)  = sgn( y,(1)  (m)y^ (m- 1)  ),  ( 6.7 ) 

where 
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y[l)  (m)  = (m)  Re [r'^  (m)]  . (6.8) 

In  order  to  use  the  modified  form  of  the  adaptive  receiver,  the  coefficients  of  the 
Lth  order  linear-predictor  must  be  calculated.  The  coefficients  of  the  linear  predictor,  a, 
are  determined  by 

a = C-1  v , (6.9) 

where 

C = B + feil/Ar0)-1I  (6.10) 

{B  lj=Rc((i-j)Tb)  (6.11) 

(6.12) 

Rc(r)  = l-{nfdtf.  (6.13) 


Figure  6.3  : Block  diagram  of  the  modified  MMSE  receiver  with  differential  detection 
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Now,  the  Lth  order  linear  prediction  of  the  channel  (amplitude  and  phase)  during  the 
current  bit  interval  can  be  formed  using  the  previous  channel  measurements  as 

(6.14) 

i=l 

where 

PS = di(m)w2]r  (m)r<1)r  (to) . (6. 15) 

The  phase  estimate  is  then  found  by  (m)  = Zf}$  (m) . 

In  order  to  adapt  the  filter  weights  in  this  chapter,  we  will  use  the  differential  least 


square  (DLS)  algorithm  [54]  as  follows: 

d'l(m)  = dl(m)yt1\m- 1)  (6.16) 

P(1)  (m  + 1)  = (1  - v)d[(m)  Re[rjj(1)  (to)]+  op(1)  (to)  (6. 17) 

R(1)  (m  + 1)  = (1  - v)  Re[r;4(1)  (m)]Re[r1,1(1)  (m)]T  + uR(1)  (m)  (6. 1 8) 

wg  (m  + 1)  = R(1)-1  (to  + l)p(1)  (to  + 1) , (6. 19) 


where  p(1)(m  + l)  and  R(1)(m  + 1)  are  estimates  of  Refr^0  (m)]}  and 

^Refr'g  (m)]Re[r1/1(1)  (m)f  } , respectively,  and  v is  a constant  close  to  one,  chosen  to 
discount  past  data.  We  use  u = 0.998 . 

6.2  The  modified  MMSE  receiver  with  multiple  antenna  sensors 

Since  the  multiple  antenna  sensors  providing  spatial  diversity  can  be  used  to 
mitigate  the  effects  of  fading  of  the  desired  signal  and  suppress  the  interfering  signals,  the 
multiple  MMSE  receivers  exploiting  spatial  diversity  may  be  expected  to  reduce  the 
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problem  of  locking  in  phase  to  any  of  the  interfering  signals.  The  possible  block  diagram 
of  the  MMSE  detector  with  multiple  antennas  with  equal  gain  differential  combining  on  a 
single  path  flat  fading  channel  is  shown  in  Figure  6.4. 


Figure  6.4  : Block  diagram  of  multiple  MMSE  receivers  with  multiple  antenna  array 

on  a single  path  flat  fading  channel 
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It  is  shown  in  [55]  that  the  conventional  MMSE  receiver  is  sensitive  to  the  fading 
rate  on  a frequency  selective  fading  channel  Very  slow  fading  rates  are  required  for  the 
MMSE  receivers  to  effectively  combine  the  multiple  paths  of  fading  channels.  In  [78],  it 
is  pointed  out  that  the  MMSE  receiver  needs  adaptive  algorithms  which  can  adequately 
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Figure  6.5  : Block  diagram  of  multiple  MMSE  receivers  using  spatial  diversity  via  two 
antennas  with  equal  gain  differential  combining  for  two  paths  frequency  selective  fading 

channels  (Structure  I) 
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track  the  time-variations  induced  by  frequency  selective  fading  channels.  One  possible 
block  diagram  to  combine  two  paths  of  fading  for  two  antennas  with  precompensated 
phase  process  is  shown  in  Figure  6.5,  which  is  called  multiple  MMSE  receivers  with 
equal  gain  differential  combining  on  two  paths  frequency  selective  fading  channels.  We 
can  easily  extend  the  structure  to  the  case  of  multipath  and/or  multiple  antenna  array. 
This  structure  is  referred  to  as  Structure  I. 

For  the  lower  bound  of  the  simulated  BER  performance  of  Structure  I in  Figure 
6.5,  the  solution  of  the  MMSE  Weiner-Hopf  equation  is  used.  In  other  words,  given 
perfect  knowledge  of  the  fading  parameters,  the  phase  variations  for  each  path  of  the 
desired  user’s  signal  are  removed  for  the  desired  portion  of  the  received  signal  by 


That  is,  we  assume  the  knowledge  of  the  channel  phase  dff  (m) . Then,  let  the 
precompensated  received  signal  for  the  rth  path  be  r^p)(m)  whose  real  part  is  the  input 
signal  of  the  modified  MMSE  receiver. 

Then,  every  tap  weight  vector  of  multiple  MMSE  receivers  with  perfect  phase  recovery  is 
determined  by  the  solution  of  the  Weiner-Hopf  equation 


forming 


(6.20) 


(6.21) 


where 


107 


2 „(rr 


= «i(.r)2c1cf  + X (au  + TU'  - T i,  VI.  + *iy  ~ tu ) 

r'= 1 4 

rVr 

+ (au(^;  + *V  -Tl  r )a^2(ATTc  +ru-  -Tlir)) 

K 2 y<-P\2 

•f££~T"Vat,i(^7’c  +T*y  -TuVwlw.  + W -t,,) 

t=2  r'= 1 4 

+ a*,2(^7;  +T*y  -Tlir  VmOvT*  + **/  -TU))+<72I 
Pu  =fi{d1'Re[r1,(/)]}  =o#)c1. 


(6.22) 


(6.23) 


where  =a($JPk/Pl  and  the  time  variable  m is  omitted. 


We  might  consider  other  lower  bounds  which  uses  different  structures,  based  on 
the  same  number  of  multiple  antennas  as  in  Figure  6.5.  The  structures  as  shown  in  Figure 
6.6  and  6.7  use  a single  (ideal)  filter  per  antenna  element  and  an  ideal  diversity  combiner 
followed  by  a single  (ideal)  filter  for  each  path  of  the  desired  user’s  signal,  respectively. 
These  structures  are  referred  to  as  Structure  n and  HI,  respectively.  Assuming  the  perfect 
knowledge  of  the  fading  parameters,  the  comparison  of  three  different  structures  by 
simulation  results  using  the  MMSE  solution  gives  us  some  indication  of  how  the  different 
structures  effect  the  performance  of  our  system. 

In  Figure  6.6,  the  tap  weight  vector  is  determined  by  the  solution  of  the  Weiner- 
Hopf  equation 

w^RW'pW  (6.24) 
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Since  the  diversity  combining  is  known  to  mitigate  the  effects  of  the  desired 
signal  due  to  fading  and  suppress  the  interfering  signals,  the  MMSE  detector  employing 
the  diversity  combiner  may  be  expected  to  reduce  the  problem  of  locking  in  phase  to  any 
of  the  interfering  signals.  In  Figure  6.7,  the  tap  weight  vector  of  the  diversity  combiner  is 

chosen  to  minimize  the  MSE,  / = £^|  d[  (m)cl  - wfp)  ( m ) Re  [r1'(/)  (m)]  | 2 j . That  is,  the  tap 

weight  of  the  diversity  combiner  is  determined  by  the  solution  of  the  Weiner-Hopf 
equation 


where 


w(p)  =R(p)~l  n(p) 
wl,r  A1  ,r  Pl,r  » 


R‘;>=£{Re[r1'j«]rRe[r';'>]} 

=<*!?2  cf«.  +j£al'/(“l(HTc  +r,,.  -ru)a  u()VTe  + r,,.  -t„) 

4 r'=l 

rVr 

+ au (fiTc  +hy  +hy  -*0) 

+IIl^/2(a‘.i(Af7’c  +*ty  -*i,r)*kANTc  +rky  -tu) 

*= 2 r'= 1 H 

+ *TkANTc  +rky  -h.rWAM'c  +zky  -tu))+<t2N 
Pu  = E{d[c\  Re[r'^)]}  = Na(xpr\ 


(6.27) 


(6.28) 


(6.29) 


109 


V 


r-  1 


1.1 


Chip 
matched 
filter  & 
store  N 
samples 


Phase  estimator 


— ►Trco  — 


L T 


1.2 


Chip 
matched 
filter  & 
store  N 
samples 


( 

>> 

Adaptive 

Differential 

Phase  estimator 

Filter 

detector 

►00— ► Re( ) 

V. 


1.1 


Chip 
matched 
filter  & 
store  N 
samples 


Phase  estimator 


L T 


1.2 


Chip 
matched 
filter  & 
store  N 
samples 


Mg>— ► 

Re() 

c: 

Phase  estimator 

A i 

Re() 

Differential 

detector 


dx  (m) , data  out 


sgn() 


Figure  6.6  : Block  diagram  of  multiple  MMSE  receivers  using  spatial  diversity  via  two 
antennas  with  equal  gain  differential  combining  for  two  paths  frequency  selective  fading 

channels  (Structure  II) 
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Figure  6.7  : Block  diagram  of  multiple  MMSE  receivers  based  on  diversity  combiner  via 
two  antennas  for  two  paths  frequency  selective  fading  channels  (Structure  HI) 


Ill 


Then,  the  tap  weight  vector  of  the  MMSE  filter  following  the  diversity  combiner  is  also 
determined  by  the  solution  of  the  Weiner-Hopf  equation 

W2,r  =R2!rP2,r  » (6-30) 
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6.3  Numerical  Results 

The  transmitted  powers  of  all  users  were  taken  to  be  equal  (i.e.,  Pk=Pi).  For  each 
data  point,  the  MMSE  filter  was  trained  for  a period  of  100  b followed  by  1900  b of  real 
data  in  which  was  updated  in  decision-directed  mode.  In  these  simulations,  the  data  rate 
is  assumed  to  be  9600  bps.  This  process  was  then  repeated  for  50  different  trials.  We 
assume  a two  path  channel  for  frequency  selective  fading.  The  paths  from  each  user  are 
assumed  to  have  the  same  mean  power.  The  third-order  predictor  for  the  channel 
estimation  is  used.  We  assume,  for  simplicity,  that  the  T!  2 is  4TC. 
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Figure  6.8  and  6.9  show  the  BER  of  the  MMSE  receiver  with  fd  =100//z  and 
Eb/N0  =15 dB  where  Eb  is  the  sum  of  the  energy  from  all  paths  of  the  desired  user’s 
signal  in  a single  path  flat  and  two-paths  frequency  selective  Rayleigh  fading  channels, 
respectively,  as  a function  of  the  number  of  users  for  three  different  number  of  antennas. 
Simulation  results  of  the  conventional  MMSE  receivers  using  (complex)  DLS  alone  used 
in  [54]  without  the  phase  precompensation  process  of  linear  prediction  phase  estimator 
(LPPE)  and  their  lower  bounds  using  the  (complex)  MMSE  solution  are  included  in 
Figure  6.8  and  6.9.  It  is  shown  that  the  difference  in  BER  performance  between  (real) 
DLS  with  LPPE  and  lower  bound  using  (real)  MMSE  solution  is  not  much.  The  DLS 
with  LPPE  can  improve  the  performance  of  the  modified  MMSE  receiver  as  compared  to 
that  of  the  conventional  MMSE  receiver  when  using  DLS  without  LPPE,  especially  in  the 
presence  of  frequency  selective  fading  channels.  As  the  number  of  antennas  inceases,  the 
performance  can  be  significantly  improved.  The  spatial  diversity  via  multiple  antennas 
can  improve  the  BER  performance  of  the  MMSE  receiver  with  a single  antenna.  It  is  also 
shown  that  the  number  of  users  the  MMSE  receiver  can  support  is  increased  by  the  use  of 
spatial  diversity  via  multiple  antenna  arrays.  In  the  simulation,  the  BERs  are  0 when  K=1 
and  K=5  for  P=4  on  the  frequency  selective  fading  channel.  That’s  why  the  zero  points 
are  not  plotted  in  the  log  scale  of  Figure  6.9. 

In  Figure  6.10,  the  BER  performance  using  DLS  with  LPPE  is  plotted  as  a function 
of  the  fading  rate  with  K=10  and  Eb/N0  =10 dB  on  a two-path  frequency  selective  fading 
channel.  The  speed  of  the  fading  seems  to  have  very  little  effect  on  the  performance  of 
the  modified  MMSE  receiver  using  DLS  with  LPPE. 
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The  BER  performance  of  the  modified  MMSE  receivers  using  DLS  with  LPPE  is 
plotted  as  a function  of  the  number  of  antennas  in  Figure  6.1 1 with  K=15,  Eb/N0  =10 dB 
and  fd  =100 Hz  for  a two-path  frequency  selective  fading  channel.  It  is  seen  that  as  the 
number  of  antenna  sensors  increases,  the  BER  performance  can  be  significantly 
improved. 

In  Figure  6.12,  the  performances  of  three  different  structures  are  compared  on  two- 
paths  frequency  selective  fading  channels  with  fd  =100 Hz  and  Eb/N0  =15 dB  by  using 
perfect  knowledge  of  the  fading  process.  The  lower  bounds  of  the  BER  are  0 when  K=1 
for  P=3  and  the  lower  bound  for  Structure  I is  0 when  K=1  for  P=3. 

6.4  Conclusions 

In  this  chapter,  we  have  examined  how  the  spatial  diversity  via  multiple  antenna 
arrays  can  be  used  for  MMSE  detection.  It  is  shown  that  the  use  of  spatial  diversity  via 
multiple  antenna  sensors  can  result  in  substantial  improvement  in  BER  and  capacity 
performance  of  the  single  antenna  based  MMSE  receiver.  By  using  the  DLS  with  LPPE 
to  adapt  the  tap  weight  vector  of  the  modified  MMSE  receiver,  the  BER  performance  can 
be  improved  as  compared  to  that  of  the  conventional  MMSE  receiver  using  DLS  alone, 
especially  on  a frequency  selective  fading  channel  and  also  the  difference  in  BER 
performance  between  the  modified  MMSE  receiver  using  DLS  with  LPPE  and  lower 
bound  using  MMSE  Weiner-Hopf  solution  is  not  much.  The  DLS  algorithm  with  LPPE  is 
insensitive  to  the  fading  rate.  The  use  of  different  structures  based  on  the  same  number  of 
antenna  sensors  does  not  much  affect  the  BER  performance. 
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Figure  6.8  : BER  on  a flat  Rayleigh  fading  channel 
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Figure  6.9  : BER  on  two-paths  frequency  selective  Rayleigh  fading  channel 
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Figure  6.10  : BER  on  two-paths  frequency  selective  Rayleigh  fading  channel 
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Figure  6.1 1 : BER  on  two-paths  frequency  selective  Rayleigh  fading  channel 
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Figure  6.12  : BER  on  two-paths  frequency  selective  Rayleigh  fading  channel 


CHAPTER  7 

SUMMARY  AND  FUTURE  WORK 


In  this  final  chapter,  the  key  contributions  of  this  dissertation  will  be  summarized. 
Some  areas  of  future  research  will  also  be  suggested. 

7.1  Overview  of  Dissertation  — Summary 

This  dissertation  investigates  how  antenna  arrays  can  be  used  for  parameter 
estimation  and  detection  in  DS-CDMA  communication  systems.  An  adaptive  antenna 
array  and  diversity  combiner  based  timing  estimation  is  proposed  which  makes  use  of  the 
multiple  antenna  array’s  capabilities  to  discriminate  between  spatially  separated  users  and 
provide  spatial  diversity.  The  capacity  of  a DS-CDMA  communication  system  using  the 
conventional  receiver  is  limited  by  the  timing  estimation  problem  and  some  multiuser 
receivers  require  accurate  synchronization.  Even  if  the  adaptive  antenna  array  and 
diversity  combiner  have  enjoyed  much  attention  in  order  to  increase  the  performance  and 
capacity  of  the  DS-CDMA  system,  they  have  received  less  attention  for  the  code  timing 
estimation  problem.  Multiple  MMSE  detections  scheme  exploiting  spatial  diversity  is 
also  proposed  on  time- varying  fading  channels. 

A general  spatial  channel  model  including  spatial  information  such  as  DO  As  and 
spatial  diversity  was  derived  in  Chapter  2.  We  then  showed  how  the  general  model  can 
reduce  to  several  specific  cases.  This  is  necessary  in  order  to  evaluate  the  beamforming 
and  diversity  combining  technique  proposed  later  on  in  this  dissertation.  If  there  is  no 
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spatial  information  like  DO  As  in  received  signals,  the  system  model  for  beamforming 
becomes  that  of  the  diversity  combiner. 

In  Chapter  3,  we  described  two  beamforming  techniques  for  tap  weight  vector 
estimation  of  the  adaptive  antenna  array  and  only  one  of  them  was  used  for  the  diversity 
combiner.  The  DMI  algorithm  was  considered  to  calculate  the  tap  weight  vector  for  the 
training  sequence  based  beamforming  or  diversity  combining  due  to  its  fast  convergence. 

Chapter  4 derived  the  Cramer-Rao  Bound  for  an  AWGN  channel  and  fading 
channels.  Two  CRBs  for  an  AWGN  channel  were  presented  and  compared  in  order  to 
investigate  how  much  information  about  the  signal’s  timing  is  lost  in  the  process  of 
beamforming.  Our  approach  for  estimating  a propagation  delay  was  shown  to  be  a good 
propagation  delay  estimation  method  under  a nearly  perfect  power  control  condition. 
Even  if  our  focus  is  on  the  code  timing  estimation  problem,  the  CRBs  in  Chapter  4 can  be 
used  as  a optimality  criterion  for  the  other  parameter  estimation. 

The  code  timing  estimation  problem  for  an  AWGN  channel  and  fading  channels 
was  dealt  with  in  Chapter  5.  The  LSML  timing  estimator  which  follows  the  adaptive 
antenna  array  and  diversity  combiner  was  proposed  for  estimating  the  propagation  time 
delay  of  a desired  signal  This  chapter  was  divided  into  three  parts  on  the  basis  of  which 
channel  model  was  assumed.  Firstly,  it  was  found  for  an  AWGN  channel  that  the 
adaptive  beamformer  could  improve  the  performance  of  the  LSML  timing  estimator  as 
compared  to  a single  antenna  based  LSML  timing  estimator.  Our  approach  could  be  a 
good  timing  estimator  even  when  a large  number  of  users  are  present  in  a perfect  power 
control  scenario.  The  CRB  obtained  in  Chapter  4 was  also  used  as  a performance  bound 
for  the  accuracy  of  the  timing  estimator.  Secondly,  the  diversity  combiner  based  LSML 
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timing  estimator  was  applied  to  a single-path  time-varying  flat-fading  channel  and 
compared  with  the  MASE  algorithm.  It  was  shown  that  it  improved  the  performance  of 
the  LSML  timing  estimator  with  less  computational  complexity  relative  to  that  of  the 
MASE  algorithm.  Thirdly,  the  efficient  algorithm  of  an  approximate  ML  approach  for 
joint  estimation  of  multipath  timings  of  a desired  user’s  signal  on  a frequency  selective 
fading  channel  was  developed  and  the  MUSIC  algorithm  exploiting  spatial  diversity  was 
also  proposed.  Their  performances  were  compared  via  simulation  results. 

Chapter  6 has  studied  the  MMSE  detection  with  multiple  antenna  sensors.  The 
spatial  diversity  via  multiple  antenna  sensors  can  improve  the  performance  of  the  MMSE 
receiver  on  flat  and  frequency  selective  fading  channels. 

7.2  Contributions 

The  main  contributions  made  in  this  dissertation  were  : 

• Developing  a general  system  model  which  includes  spatial  information  based  on 

multiple  antenna  sensors 

• Studying  beamforming  or  diversity  combining  technique  for  DS-CDMA  signals 

• Derivations  of  the  Cramer-Rao  bound  for  both  AWGN  and  fading  channels  when 

using  multiple  antenna  array 

• Comparing  the  post-combining  CRB  with  the  pre-combining  CRB  on  AWGN 

channel 

• Developing  timing  estimators  for  both  AWGN  and  fading  channels 

• Developing  multiple  MMSE  receivers  with  multiple  antenna  sensors 
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7.3  Future  Work 

Some  topics  of  future  research  as  extensions  to  the  problems  studied  in  this 
dissertation  will  be  suggested  in  this  section. 

We  have  considered  only  the  training  sequence  based  diversity  combining  on  the 
time-varying  fading  channel  because  the  DOA  based  beamforming  may  not  be  effective 
due  to  the  time- varying  phase  of  the  fading  process.  Future  work  will  therefore  include 
the  study  of  effective  algorithms  for  estimating  the  DOA  of  the  desired  user’s  signal  and 
also  effective  beamforming  algorithms  on  the  time-varying  fading  channel.  Further, 
modifications  to  both  the  training  sequence  based  beamforming  and  the  DOA  based 
beamforming  for  frequency  selective  fading  channels  may  be  required. 

Another  significant  area  of  future  work  is  the  estimation  of  the  fading  processes  of 
the  CDMA  signals.  It  is  shown  in  [78]  that  the  conventional  MMSE  receiver  will  have 
substantial  degradation  in  performance  if  it  is  not  able  to  track  the  fading  channel 
parameters.  We  have  used  linear  prediction  phase  estimator  (LPPE)  based  on  DLS 
algorithm  for  the  modified  MMSE  receiver  to  estimate  the  fading  process  of  each  path  of 
the  desired  user  at  each  antenna  on  a single  path  flat  and  two  path  frequency  selective 
fading  channels.  Even  if  the  DLS  algorithm  with  LPPE  could  properly  track  the  fading 
process  of  the  desired  user’s  signal,  the  DLS  algorithm  requires  the  inverse  operation  of 
the  autocorrelation  matrix.  Therefore,  it  is  of  interest  to  study  more  efficient  adaptive 
algorithms  in  terms  of  computational  complexity  while  achieving  the  MMSE  solution. 

Our  work  in  this  dissertation  focuses  on  using  multiple  antennas  at  the  receiver  to 
improve  the  quality  of  the  received  signal.  Recently,  the  transmit  diversity  schemes  using 
multiple  transmit  antennas  have  been  shown  to  reduce  the  effect  of  fading  at  the  remote 
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units.  Therefore,  if  we  apply  our  techniques  developed  in  this  dissertation  to  system 
model  with  transmit  diversity  as  well  as  receiver  diversity,  the  system  performance  could 
be  improved.  It  is  of  interest  to  study  the  performance  of  system  with  transmit  diversity. 

Another  area  that  is  receiving  attention  in  the  literature  is  that  of  multicarrier  DS- 
CDMA  systems.  In  a multicarrier  DS-CDMA  system  a data  sequence  multiplied  by  a 
spreading  sequence  modulates  multiple  carriers  rather  than  a single  carrier.  It  has  been 
shown  that  a multicarrier  DS-CDMA  system  is  robust  to  both  narrowband  interference 
and  multipath  fading.  However,  it  requires  the  synchronization  between  the  spreading 
pseudonoise  code  of  the  desired  user’s  signal  and  the  despreading  code  locally  generated 
in  the  receiver.  Therefore,  it  is  of  interest  to  study  the  performance  of  synchronization 
and  detection  of  our  approaches  developed  in  this  dissertation  for  multicarrier  DS-CDMA 
communication  systems. 


APPENDIX  A 

DERIVATION  OF  THE  CRAMER-RAO  BOUND 

The  derivation  of  the  CRB  for  the  AWGN  channel,  in  Section  A.1,  and  for  the 
fading  channels,  in  Section  A.2,  when  using  the  multiple  antenna  sensors  is  an  extension 
of  the  CRB  based  on  a single  antenna  sensor  presented  in  [47]  and  follows  closely  the 
procedure  presented  in  [48] 


A.1  CRB  for  the  AWGN  Channel 


We  start  by  repeating  some  key  equation  from  subsection  4.1.2  of  Chapter  4.  The 
received  vector  for  the  mth  data  interval  from  P antenna  sensors  is 


r(m)  = B(r)E(y,  d)z(m)  + n(/n)  , 


(A.1) 


where  n(m)  is  a complex  Gaussian  random  vector  with  zero  mean  and  second  moment 


matrices 


n(0n  "(;')]  =cr2I  NP8K(i-j) 
n(0nr  (;')]=  0. 


(A.2) 


Conditioned  on  'P  and  z(m) , the  received  vector  r(m)  will  be  complex  Gaussian  vector 


with  mean 


£’[r(m)|'P,z(m)]=  BEz(m) 


(A.3) 


and  second  moment  matrices 


124 


125 


4(r( 0 - BEz(OX?(;)  - BEz(;)f  \¥,  z(m)]  =o2lNP8K(i-  j ) 
E[(?(t)  - BEz(i)X?(;)  - BEz(;)  f fP,  z(m)]=  0 . 


(A.  4) 


The  probability  density  function  for  r(m)  conditioned  on  4*  and  z (m)  is  called  the 
likelihood  function  for  r(m)  and  is  denoted  by  L(r(m)) . 

The  likelihood  function  of  r (m)  can  be  written  as 


L(f(m))%-de,^/„rP('^11  ?(m)-  BEz<m)i! 

— tII  ?("0-BEz(m)||2 


(A.5) 


Since  r(i')  and  r(j)  are  uncorrelated  complex  Gaussian  random  vectors  for  i * j , the 
likelihood  function  of  r = [ Fr  (1)  rr  (2)  ?r  (Af )]  may  be  written 

L(r)  = n^(m)). 


as 


(A.  6) 


m=l 


By  taking  the  natural  logarithm  of  L(r(m))  and  ignoring  terms  that  are  independent  of 
4* , the  log-likelihood  function  is 

1 M 

In  L(r)  = -MNP In  cr2 r£||  f(m)-BEz(m)|  . (A.7) 


m- 1 


A.  1.1  Differentiating  the  Log-Likelihood  Function 

We  start  now  computing  the  gradient  of  the  HM,  9ln  L(r)/34/ , as 


d In  L(r) 


d In  l(?) 


do2 

d In  L(r) 


dy 

<9lnZ-(r) 


dd 

d In  L(r) 


. dr 


(A.  8) 
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The  Fisher  information  matrix  (4.24)  can  be  expressed  as  4 x 4 block  matrix  where  each 
block  is  the  expected  value  of  the  outer  product  of  two  gradients. 

In  order  to  compute  the  gradients,  we  will  use  the  following  results. 

RESULT  1 : Let  x e cn  be  a function  of  0 = ^ ...  jr  e Rm , and  suppose 

A e C"*"  is  a constant  matrix.  Then 


— x"Ax 
50 


Ax  + 


* AV. 


(A.9) 


where 


dx  _ dx  dx 
50  ~~  |_5^~  ~d02 


eC" 


(A.  10) 


RESULT  2 : Let  x e C"  be  a function  of  0 e Rm , and  suppose  a e C"  is  a constant 


vector.  Then 


(A.  11) 


(A.  12) 


RESULT  3 : Let  x e Cnx*  be  a function  of  q e r , and  suppose  a e C*  is  a constant 
vector.  Then 


— a"X"Xa 


dO 


a 


H 


aC\H 


Xa+a^X" 


faC\ 
^ ) 


a- 


(A.  13) 
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i)  Computing  d In  L(r)/3cr2  : 

We  start  by  differentiating  the  log-likelihood  function  with  respect  to  a 2 

d In  L(?)  MNP  1 A.*,  x 

a 2 = — + — Xn  (»*)n(m) . 

do  O ~ 


m=l 


ii)  Computing  3 In  L(r)/3y  : 

Straightforward  differentiation  with  respect  to  yt  yields 


where 


Now, 


= ' — tS4-  (fo»)  ~ BEz(m)f  (r(m)  - BEz(m)) 
o *Yk 


U m=l 

1 W 

=^rS2Re 


( ze  \H 

-B— z(m) 

dyk 


n(m)  + nH  (m) 


( 


-B-^-z(m) 
&Yk 


( jp 

B— -z(m) 


ii(m) 


= ^rX2Rejz,,(m)^-BHn(m)|, 


m=l 


3EW 

-t— 0.;)= 


exp(-;0t)  if  i = j -2k -l 
exp (~jdk ) i/  i = j = 2k 
0 otherwise 


dEH 

= exp(-./0t)[z2t_1(m)  ^2*('”)] 

or* 


K# 

')2H 

bff 
2Jt  . 


n(m). 


(A.  14) 


(A.  15) 


(A.  16) 


(A.  17) 


and  the  gradient  with  respect  to  y can  be  expressed  as 
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2i2iS=-LE2Re 

dr 


expC-;0,  )[zx  (m)  z2(m )] 


a n(m) 


exp(-;0jf)[z2jr_1(ro)  z2jc(m)] 
= ~V  £ 2Re{z*  (m)E*  B*n(m)}, 


>2*-l 
Lb"x  j 


n(m) 


(A.  18) 


where  Er  e C2*x2*  is  defined  as 


E = diag(exp(jdl),exp(jdl ), ,exp(jdK),exp(jdK))  (A.19) 


and  Z(m)  e {- 1, 0,  l}2*x*  is  defined  as 


(z  (m)Xi,k)  = 


z2k- i(m)  j/  j = 2/r  - 1 
z2*  (m)  j/  i = 2k 
0 otherwise 


(A.  20) 


iii)  Computing  d In  L(f)/dd  : 

If  we  can  follow  the  same  steps  as  the  previous  one  (ii),  we  can  get  the  following 

form 


- = -Li2Re{z"  (m)E»“ B“ii(m)}, 

cW  CT  m=1 

where  E0  e c2**2*  is  defined  as 

Ee  — diag(j^ , j^i , j%2 , j%2, * j^K  * J^k  ) = 7®  • 

iv)  Computing  3lnL(r)/5r  : 

Straightforward  differentiation  yields 


(A.21) 


(A.  22) 


^ In  f,(r ) _ 1 ri 


-nH  (m)^-Ez(m)-zH  (m)EH  n(m) 
dzk  drk 


= ^-X2Re|Z,/(m)Eff  n(m)|. 


m= 1 


(A.23) 
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Recall  the  expression  for  the  common  elements  of  the  columns  of  the  components  of  B 

(A.  24) 


au  = 


a*.2  = 


8k 

f 5 A 1 

-^P(pt+U)  + 

1 ~~r  P(pv,l)  c 

Tc 

1 

. <j 

8k 

f s A 

-fP(pt+l,-l)  + 

l~P(P*.-Dc 

_ c 

(A.25) 


As  seen  from  these  equations  the  derivative  of  B with  respect  to  zk , or  equivalently  with 
respect  to  8k , might  not  exist  when  8k  = 0 . The  reason  for  this  is  that  ak  2 lies  on  a 
straight  line  between  P (pk  + 1,  - l)ct  and  P(pk , - l)ct  when  zk  e [pk  Tc,(pk  +])TC ) and 
between  P^*  + 2,  - l)ck  and  P (pk  + 1,  - l)c*  when  zk  e [(pk  + 1 )TC  ,(pk  + 2 )TC ).  Now  if 
P (pk , - l)ct , P(pk  + 1,  - l)ct  and  P(pk  + 2,  - l)c*  are  not  collinear,  then  the  derivative 
will  not  exist  for  zk  = (pk  + 1 )TC . The  CRB  is  not  valid  if  the  derivative  does  not  exist, 
(see,  for  instance,  Van  Trees  [58])  In  order  to  proceed,  we  therefore  make  the  assumption 
only  that  8k  * 0 for  Jt  = 1, 2,  •••,  K . Since  the  assumption  only  implies  that  we  are  not 
sampling  in  synchronization  with  any  user,  it  is  a rather  mild  restriction. 

Under  the  assumption  that  8k  * 0 for  k = 1,2, ,K , 


dB 

dzk 


0 


o 


db 


2k 


dzk  dzk 


0 


0 


(A.26) 


Now,  if  define 


g 


db 


2k— 1 


2k— l 


*k  ’ 


_ db2* 
gu  ’ 


(A.27) 


then. 


dRH 

iH  (m)Ew  — — ii(m)  = [z2k.x  (m)  z2k  (m)] 
Wk 


aH 

b2k-\ 

zH 

. 62 * . 


n(m). 


(A.  28) 
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and  the  gradient  with  respect  to  x can  be  expressed  as 


^ln_Lg)  = 1 f 
* 


^rUi(»j)  z2(wo]  8lff  “(«> 

Lg"  j 


ljrb2x-i(w)  z2*("*)] 
= £ 2 Refc  " (m)E  " G " ■<»>}• 


g£c-i 

Lgfx  J 


n(m) 


(A.29) 


m=l 


where  G e C^2*  is  defined  as 


G = (g,  g2  g2Jf]. 


(A.  30) 


A.  1.2  Forming  the  FLsher  Information  Matrix 

In  order  to  compute  the  expected  value  of  the  outer  product  of  all  the  gradients 
needed  to  form  the  FIM,  we  need  the  following  results. 

RESULT  4 : 


E 


[nH  (m)n(m)fiff  (fc)n(Jfc)]=  • 


(NPfa* 
NP(NP  + l)o4 


if  m*k 
if  m = k 


RESULT  5 : For  all  m,  k. 


£[ntf(m)n(m)nr(jfc)]  =0  . 
RESULT  6 : Suppose  x,  y € C* , then 

Re{x}Re{yr  }=  ^ [Re{xyr  }+  Rejxy  w ]] 
Im{x}lm{yr  }=  [Re{xyr  }-  Re{xyw  J 
Re{x}lmjy  r }=  i [lm{xyr  }-  Im{xyw  ]j . 


(A.31) 


(A.  32) 


(A.33) 
(A.  34) 


(A.35) 
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The  (1,1)  block  of  the  Fisher  information  matrix  is 

2MNP  —// , _ 1 _ h _ x—  // 

— — 2,n  (m)n(m)  + — 2^Xn  (/n)n("J)n  (*)“(*) 

^ m=l  ^ m=l  4=1 

= (AfA^) — lM0_MNPa2  +-L[mnp(NP  + 1)cj* +{]U2  - m\nP)2o 4] 

a a6  o 


(A.36) 


MNP 
■ a4  ' 

We  note,  by  using  Result  5,  that  d In  L(r)/9<r2  is  uncorrelated  with  all  other  gradients. 
Thus,  the  (1,2),  (1,3)  and  (1,4)  blocks  of  the  information  matrix  will  be  zero. 

To  find  the  remaining  blocks  of  the  information  matrix,  we  need  to  compute  the 
expected  values  of  all  possible  outer  products  of  the  gradient  (A.  18),  (A.21)  and  (A.29). 
The  amount  of  algebra  required  can  be  reduced  by  noting  that  the  gradients  are  all  of  the 
form 


4-ZRe^H(m)n(m)}, 


(A.  37) 


m= 1 


where  F(m)  is  some  NPxK  complex  matrix.  The  general  expression  for  a block  in  the 
information  matrix  is  then 


E X Re(F"  (m)n(m)}T ^ Re{F2"  (A:)n(*)}l 


<7  m=l 
M U 


{<7  *=l 


= Re{F*  ("0n(m)(F2H  (k)a(k)J  }+  Re  {f"  (m)n(m)  (f"  (k)n(k)f  }l 

"=» k=l  1 L J(A.38) 


M M 


= -4SERe  {F"  (m){<72lN8(k  -m))F2(k)} 

° m=l  i=l 
2 M r 

=— 5>e{F"(ro)F2(ro)}. 
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The  blocks  of  J,  as  found  in  (4.25)-(4.30),  follows  by  identifying  the  corresponding 
F(m)  in  (A.  18),  (A.21)  and  (A.29),  and  substituting  the  result  into  (A.38). 


A.2  CRB  for  Fading  Channels 


We  will  now  derive  the  CRB  for  fading  channels  modeled  as  in  (4.49). 


A.2.1  Differentiating  the  Log-Likelihood  Function 
i)  Computing  d In  L(r)/9<x2 ; 


We  start  by  differentiating  the  log-likelihood  function  with  respect  to  <7 2 , 

dO  GO  m=i 


(A.  39) 


ii)  Computing  d In  L( r)/dp  ( m ) ; 


Consider 


(A.  40) 


Using  Results  1,  2,  and  3,  we  can  show  that 


9lnZ,(r)_  1 
9/J(m)  a2 


f3S(m)fx„S(m)  + ['3|m)TxrS-<m) 


dy3  (m) 

( 9E(m)  ^ H 

dp(m) 


9/J(m) 


XHn(m) 


(A.41) 


where 


dZ(m) 

dfl(m) 


dZ(m)  9 Z(m) 


dpi(m)  BP  2 (m) 


9E(m) 
dpK  (m) 


eC 


2 RPxRP 


(A.  42) 
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35  (m)  _ 35  (m)  35  (m) 

3 &(m)  " 3&(1)(/n)  3&(2)(m) 


35(m) 

3^(P)(m) 


Q2RPxRkP 


35  (m)  _ 35  (m)  35  (m) 

3 ®}(m)  3®J(m) 


35  (m) 

3ft!£  (*>) 


35  (m) 


0—0 


0^0 
2(r-l)  jeroi 


3flff  (m)z2t_1  (m)  3/3ff  (m)z2t  (m) 
3ffi?(m)  dfiffOn) 


0—0  0—0 


where_^Lec!*~ 

3ft;’ (m) 


and  0 is  a 1x2/?  row  vector  of  all  zeros  and  where 


dPk'J^r^~ = aflWr  i ^ (m) + ^ (m)lZ2*-1  (m) = Z2*-1  (m) 

dptyim)  oPkr(m) 


3/?ff(m)z2*(m) 


= z2k  (m) . 


If  we  define  Zt(m)e  {-l,0,l}2*‘x*‘  and  Di(m)£R2*,M,'>  as 


(Zk(m)Xi,j)  = 


Cm)  ifi  = 2j-\ 
Zu  (m)  ifi  = 2j 
0 otherwise 


Dt  (m)  = diag{Zk  (m),  Zk  ( m ), , Zt  (m)), 


then,  we  can  easily  form  3S(m)/3/3(m)  as 


|=t”;  = diagip^m),  D2  (m), , D*  (m))= Z(m) . 

3/3  (m) 


Thus, 


4Re(z"  wx'«(4 

3/3  (m)  o’ 


(A.  43) 
(A.  44) 

(A.  45) 

(A.  46) 
(A.  47) 

(A.48) 

(A.49) 

(A.50) 

(A.51) 


iii)  Computing  3 In  L(r)/3/3(m) ; 
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Similarly  as  for  P (m) , it  can  be  shown  that 


L'giiLj.Re 

dp  (m)  cr 


3H(m)  ^ 
df}(m)  ^ 


XHn(m) 


and 


Thus, 


dE (m)_  . dE(m) 
dp(m)~J  dp (m) 


<9  In  L(r)_  2 
dp(m)  a2 


Re 


/ 3E(m) 


XHn(m) 


=-^Im{zH(m)Xwn(ro)}. 


iv)  Computing  c?  In  L(r)/3r ; 

Using  Results  2 and  3,  we  can  show  that 


d In  L(t  )_  1 

^k,r  O2  h 


ZT(m) 


r_dxy 


n'(m)  + ZH  (m) 


( ax  ' 


n(m) 


0 M 


m=l 


/ \ff 

' ax  ' 


an 


V *-r ) 


n(m) 


Consider 


ax 


0-0 

(k-\)  matrices 


dTk,r 


0 — 0 


eR 


NPx2RP 


where  0 is  a NP  x 2Rt  P zero  matrix  and 


ax. 


= diag 


d\k  d\k 


d\k 


fok.r  ’ d*k.r  ’ 


€ R 


(A.52) 


(A.53) 


(A.  54) 


(A55) 


(A.56) 


(A.57) 
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9At 

dTk, 


0...Q  0 -0 


2(r — 1)  column  s 


d*k,r  d*k.r 


e R 


Nx2Rk 


We  define  a'2r_, , a*  2re  R**1  as 


_ ^**,21-1  _ d*k,2r 

w = ^ = 3r,r 


and  rt(m)e  C2**™*  and  Ht  e Rwv2*‘/>  as 


(A.58) 


(A.59) 


(r k(m)Yi,j)  = < 


PlyWtu-  i(m) 
flky  Mix  (m) 


0 


if  i = 2 j - i 
if  i = 2 j 
otherwise 


(A.  60) 


where 


P = 


i-l 

2^ 


+ 1 and  / = j-nxRk  where  n is  the  integer  part  of  the  quotient, 


i — - and 
Rk 


H*  = diag(\r 


(A.61) 


A 


d\k 


Now,  we  can  form 


c?lnL(r)_  2 
dz  a2 


XRe{r”(m)H"n<m)}, 


where 


r(m)  = diag(Tx (m),  r2(m),  , TK(m))e  C2RFxR 

H = [H,  H2  HjeRMiP. 


(A.62) 


(A.63) 


(A.  64) 
(A.65) 


Summary 
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As  shown  above,  the  derivatives  of  the  log-likelihood  function  can  be  written  as 


da  a a m=i 

(A.66) 

dp{m)  a 

(A.67) 

dp(m)  a2  1 

(A.68) 

3ta/,(r-^TRe{r“  (•)»'■(■))' 

dx  a2  £ 

(A.69) 

A.2.2  Forming  the  Fisher  Information  Matrix 

By  following  the  similar  calculations  used  in  [47]  in  order  to  get  the  expected  value 
of  all  possible  outer  products  of  the  matrices  (A.66),  (A.67),  (A.68),  and  (A.69),  we  form 
the  blocks  of  J. 


i)  E 


d In  L(r)(  d In  L( F)  Y 


da2 


da2 


MNP 


(A.  70) 


ii)  By  using  Result  5,  it  can  be  shown  that  d In  L(r)/da2  is  uncorrelated  with  all  other 
derivatives.  Thus,  the  (1,  n)  blocks  (and,  therefore,  the  (n,  1)  blocks)  of  J will  be 
zero. 


iii)  E 


<?lnL(r)r  dlnZ^r)^ 
dp(m)  dp(m) 


= 4-Re[z"  (m)XwXZ  (p)]&r  (m  - p) 


(A.71) 


J5 

i 

^lnL(r)> 

T~ 

dp(m) 

v dP(m)  j 

-- ^Im[z*  (m)XHXZ(p)^K ( m- p ) . 
<T 


iv)  E 


(A.72) 
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dlnZ^iOfdlnZ^r)^ 
E — ^ 


V) 


dp(m) 


K dp(m) 


= -\] Re[z"  (m)XHXZ(p)^K  (m  - p) 


= E 


d In  Z,(r)f d In  L(f)^T 
dj5(m)  ^ dfi(m)  ; 


vi)  E 


<7lnL(rY  dlnZ^r)^ 
dp  (m)  dz 


= -^Re[zw(m)X*HT(m)]. 


vii)  E 


dlnL(rV  <?lnZ,(f)Y 

dp(m)  K dz  j 


= -^Im[zH(m)XHHT(m)]. 
cr 


viii)  E 


d In  L(r  ) f d In  L( r ) Y 


<9t  [ dr  , 


= ^-Re[rff(m)HtfHr(m)]. 


(A.73) 


(A.74) 


(A75) 


(A.76) 


We  are  now  ready  to  form  the  FTM  J.  The  information  matrix  is  partitioned  as 


where 


J = 


**  OO 

Ji  J 

JL  Jr 


&P 

PP 

T 

Pr 


'px 


Jm=E 


d]nL(?)f  d\nL( r)^ 
da 2 dcr2 


MNP 


J a*=£ 


JOT=£ 


t?lnL(r)f  <7lnZ,(r) 


J~=£ 


<9lnL(r)T  d In  L(r)  Y 


M 


= ^-ZRe[rw(m)H,/Hr(m)] 

O'  m=l 


(A.77) 

(A.78) 
(A.79) 
(A  80) 
(A.81) 
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and  where  and  J^T  are  defined  as 


J pp  — E 


c?lnL(r)f  <9lnL(r) 


\T 


j 0T  ~ E 


din  L,(f)(  dlnZ,(r)Y 
dp  I dr 


(A.  82) 


(A.83) 


From  (A.71),  (A.72),  and  (A.73),  we  note  that  the  derivative  of  lnL(r)  with  respect  to 
P(m)  and  p(p)  are  uncorrelated  when  m*  p.  1 pp  will  therefore  be  block  diagonal 

J 00  = diagifi 00  (1),  J pp  (2),  , J pp  (A/)),  (A. 84) 


where  the  m th  block  is 


J 00  (m)  = 2 


_2 

a 


Re{z/,(m)X*XZ(m)}  -Im{zH(/n)X"XZ(m)] 
-Imrfcff(m)X*XZ(m)}  Re{z"  (ro)XffXZ(ro)j 


Furthermore,  from  (A.74)  and  (A.75)  we  see  that 

J„=[j *0)  J*<2)  Jj,(M)]r, 


where 


JflT(m)=— 

Cj 


Re 

Im 


Zw(/n)XHHT(m) 

ZH(/n)XtfHT(m) 


In  summary,  we  can  write  the  information  matrix  as 


J = 


J oa  0 0 

J 00  J 0T 

J L J tt 


9 oo 

o 

o 


(A.85) 


(A.  86) 


(A.  87) 


(A.  88) 


A.2.3  Inverting  the  FIM 
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The  CRB  matrix  is  found  as  the  inverse  of  the  Fisher  information  matrix.  Since  we 
are  mostly  interested  in  bounding  the  performance  of  a propagation  delay  estimator,  the 
bound  on  the  covariance  matrix  of  such  an  estimator  is  obtained  as  the  CRB  matrix.  In 
other  words,  we  can  write  the  CRB  matrix  as 


J!  = 


CRB  (a2)  0 0 

0 CRB(p)  x 

0 x CRB{  T) 


(A.89) 


where  x denotes  matrices  that  are  of  no  particular  interest  to  us  and  where 


(A.90) 


for  any  unbiased  estimator  f of  the  propagation  delay  vector  r . 
By  using  the  results  in  [47],  we  can  find  the  desired  CRB  as 


crb~\t)=jti  -j^jSj* 

= ^y>e[r"(m)H"HT(m) 

o t* 

(m)HHXz(zH(m)XHXZ(m))"1ZHXHHr(m)J  (A'91) 

- 4 f. KelT"  ('")H"Pw.)HT(m)]. 
a *= i 


where  Pxz(m)  is  the  projection  matrix  onto  [range  (XZ(m))]1 , 


Pico  = -XZ(m)(zw  (m)Xff XZ(m))1  ZH  {m)XH  . 


(A.  92) 


APPENDIX  B 

DERIVATION  OF  EQUATION  (5.33) 


From  the  definition  of  au(ru) , au(r) , a12(ru)  and  ali2(r) , we  have 

au(Tu)  = [a,,i(pi,irc)  a^- ( [pi j + l]rc  ^ 

a u (r)  = [a,,  (PTC ) aw ( [p  + 1 ]tc )]^  “ ^ j 


and 


1-Mu 

a;,  =[ai,(pu7;)  aw(tu+l]re)  au(pre)  au([p+l]re)]  ^ 

0 


0 

0 

1 -M 
M 


and 

Au*Q‘^kj(Pu^)  au(tu+l]re)  a u{pTc)  au([p  + l]re)] 
where  i = 1, 2. 


And  then. 


where 


AuQ'X  +A$Q-1A;i2  =Ur(A"Au  + A*  A1>2)u  = UrHU, 


1 M14 


Mu 

0 

0 


0 

0 

1-M 
M . 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 
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Then, 


H = \H A +\H A 

11  “ A1,1A1,1  ^ a1,2"1,2 


#11 

#12 

#13 

# 

#2. 

#22 

#23 

# 

#31 

#32 

#33 

# 

#41 

#42 

#43 

#, 

A14  Q_1Ai4  + Aj^Q_1A{2 


*10  + *llMU  + *12^14  *20  + *21 /hj  + *22^uM  + *23/4 

*30  *3lMl4  *32/4l,l/4  + *33M  *40  *4lM  *42M 


where 


*10  — #11  * *11  — ~ ^#n  #21  #12  » *12  = #11  — #21  — #12  #22 

*20  — #13  ’ *21  = _ #13  #23  » *22  = #13  — #23  _ #14  #24  » *23  = ~ #13  #14 

*30  = #31  * *31  = _ #31  #32  » *32  = #31  “ #41  ~ #32  #42  » *33  = _ #31  + #41 

*40  = #33  » *41  = _ ' 27/33  ■*"  #43  ■*■  #34  * *42  = #33  - # 43  ~ #34  #44  • 

Then, 

(auQ-'Au  +a;jq-a;j)“  =^- , 

where 

jy  _ *40  + *4lM  + *42M  *20  *2l/4l,l  *22/4l,l/4  *23/4 

*30  *3l/4l.l  *32/4lj/4  *33/4  *10  +*ll/h.l  ^ *12 Ml 4 

D = rf0  +d1/iu  + d2/i  + rf3/i11/r  + d4/r21  +d5/i2  +rf6/x,4/r  + rf7/r14/r2  +rf8MuM2 

where 


^0  — *10  *40  *20  *30 

= *44*40  “*20*31  —*21*30 


(B.7) 


(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 


^2  *10  *41  *20  *33  *23  *30 


(B.18) 
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And, 


where 


^41  “^20^32  ~^21^33  “^22^30  “ ^23^31 

(B.19) 

^4  = ^12^40  “^21^31 

(B.20) 

^5  = ^10^42  — ^23^*33 

(B.21) 

^6  = ^12^41  — ^21^32  “^22^31 

(B.22) 

^7  = ^11^42  ~ ^22^33  “^23^32 

(B.23) 

^8  = ^12^42  _ ^22^32  * 

(B.24) 

AyQ'1d[,)  + Ai®Q_1d(/)  =Ur(A1^d1(")  +A1fl2d^))=Urg(/,),  (B.25) 

d^)  = Q~^d\p) , d(2p)  =Q"^d(/)  (B.26) 

g(p)=A1w1d1(p)+A1w2df  =\glp)  g[p)  g™  giP)]T.  (B.27) 


Then, 


where 


Then, 


where 


A-Q-‘d<'>  + A;"Q_1d^'  = 


'H, 


-i  A</>  = 


q[p)  +q(2p)Hu 

. q(p)  + q\p)n 


(B.28) 


q[p)  = g[p) . q[p)  = -glp) + g[p)  * = ^3P) . qT  = -g\p)  + gip)  - 


(B.29) 


(A'fo-'A;J+A'“Q-|A;iJ1(A;J»Q-'d“+A;“Q-ia<'') 


/,‘.w  + /,','V,  4 +fS,n+fg’M+fX'i? +fg’nun 


•(/>), 


-0>) 


(/>), 


•(/>),  ,2 


Kj>), 


°L/»)+/a%  +/a  V + /b  Vu/i  + /m 'M2  + /a  Vu/i2. 


-O’) , 


-O’) 


(B.30) 


/,SW  = A.®,1"  /iw  = -Ax,®;'1  + A ,„?S 


,00 


,00 


(/>)  _ 


,(p) 


(B.31) 
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Thus, 
C 


'3,2 

P 


P- 1 
P 

V 

— P= 1 


h(')-*J0^)+*U^3'> 

(B.32) 

0 = -h33q[P)  + h10q(4p) 

(B.33) 

~^hiQ\P^  ~ +huq(/ 

(B.34) 

~^hlQ2P)  +^12<73P) 

(B.35) 

-^32^2P)  ■*■  hnq*’'*  ’ 

(B.36) 

"‘A^r^uQ-'d^+A^Q- 

-'d<'>)} 

2 + "fi'VuM  + «7P)MuM2  + n^nln2 ) 

(B.37) 

" /OO 
J 20 

(B.38) 

(B.39) 

■9(/''/ir) 

(B.40) 

f(P)  i _(j)*  f(P) 
J 23  ^ Qa  J 21 

(B.41) 

“ f (P) 
J2A 

(B.42) 

7£rt 

(B.43) 

D 


where 


(B.44) 
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Hence, 

C\1  (^1,1  > 7 

i>r 


„(/>)  *<J>  , „(/■)*  f <*>  . ^OO*  f 00 

”7  “<7l  Jl5  + ?2  /l4  +<?4  /23 


nsP)  =4iP)*  fisr)  +Q4P)t  fx*  * 


(B.45) 

(B.46) 


s_  AT(^u,/j) 

d(m.i>h) 

+nlP)Hul  +n(2p)n+nip)nuin+ni/)nl1  +n^V  + n(6p)nl1n+n(7p)^lln2  +n(sp)nhn2) 


d0  +d1nll+d2n+diHiiH  + d4nll+d5n2  +d6vlln  + d7nuln2  +dsnl1n2 


(B.47) 
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